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Introduction 


In Unit 4 you met Newton's second law and the law of terrestrial gravitation. If we 
consider an object falling near the earth’s surface, and ignore forces other than 
gravity acting on the object, then these laws lead to a model of its motion in 
which the object's acceleration is constant. If the height of such an object above 
ground level is x at time ¢ then the mathematical statement of this model is 


c= 

a, = 
where vp = & is the velocity and g is the gravitational acceleration. This equation 
can also be written 

@x 

ge 78 (1) 

aes ee re (a é 

where qe san abbreviation for the second derivative alan of x with respect to 


t, Because Equation (1) contains a second-order derivative it is called a second- 
order differential equation, and as such is an example of the sort of equations that 
will be discussed in this unit. 

In Unit 4 we saw that Equation (1) can be solved rather simply, by direct 
integration. We have 


Ra ae 
a \a) ~~® 


Integrating with respect to t gives 
dx 
a a fear +A 
=—git+A 
where A is a constant. Integrating a second time gives 
x= —dgt?+Ar+B (2) 


where B is another constant. Equation (2) is the general solution of the second- 
order differential Equation (1). The constants A and B in the general solution may 
take any value—they are arbitrary constants, 


This example illustrates some general features of second-order differential 
equations. The appearance of two arbitrary constants in the general solution is 
typical of second-order differential equations. In Unit 2 you saw that one condition 
was needed to pick out a particular solution of a first-order differential equation. 
For a second-order equation, two conditions are required to pick out a particular 
solution. These conditions set values on the two arbitrary constants. In the 
example, we can find values for the constants A and B if we know the position x 


and velocity & at time t = 0. Equation (2) can then be used to predict the 
position of the object at all future times 1. 


The example also illustrates very well why second-order differential equations are 
important. They arise commonly in modelling. In particular, they arise frequently 
in dynamics, because dynamical models often involve the use of Newton's second 
law. Mathematical statements which are derived using this law involve 
acceleration, the second time-derivative of position. However, our concern in this 
unit will be with the mathematics of solving second-order differential equations. 
Elsewhere in the course you will see many examples of their use in modelling. 
We have seen that the general solution of Equation (1) can be found by direct 
integration. However, this method only works for equations of the form 

a@ 


ad 
ae =f 
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(where the integrals of f(t) and froa can be found explicitly), Most of this unit is 


devoted to describing a method of finding analytical solutions of a wider class of 

second-order differential equations. The technique does not depend on integration 
(it is easier!). Rather it is a matter of recognizing that the equation is of a certain 

type, for which a set routine then provides a solution. 


The second-order differential equations with which we are mainly concerned are 
linear constant-coefficient equations. (I will define these terms shortly.) An example 
is 
ay 

430 -y= 4 sin 2. (3) 
Equations of this type have the advantage of being relatively straightforward to 
solve analytically. 
Equations which are not of this special type may not be so cooperative. In such 


cases numerical methods may have to be used. The basis of one such method is 
described in Section 5. 


Terminology and notation 


You will need to become familiar with certain jargon associated with second-order 
differential equations. (You may recognize some of it from Units J and 2.) We 
shall be exclusively concerned with linear equations. 


A second-order differential equation is linear if it can be brought 


to the form 
dy dy 
Plxdaa + ax) Ge + r(x)y =f(x) (4) 


where p, g, r and f are known, continuous, functions, and p is not 
the zero function. 


I shall make a habit of writing linear equations in the form of Equation (4), with 
all the terms involving the dependent variable (in this case y) on the left, so that 
the right-hand side is a function of the independent variable (in this case x) only. 


The equations we will be particularly concerned with are constant-coefficient 
equations, that is, equations for which p, q and r are constant functions. Equation 
(3) above is an example of a linear, constant coefficient, equation. 


We concern ourselves here only with linear constant-coefficient equations because, 
as I remarked earlier, these are relatively easy to solve. However, even when the 
coefficients are not constant the idea of a linear differential equation is still 
important because there are a number of general theorems about differential 
equations which apply to linear equations. 


There is one last piece of terminology which we shall need when we study methods 
of solution. An agp which can be brought to the form 


nix) © y+ ae & + rx)y=0 


(that is, —_ (4) with the zero function on the right-hand side) is a 
homogeneous linear differential equation. If the right-hand side of Equation (4) is 
non-zero, then we have an inhomogeneous equation. 


Exercise 1 
(i) | Write down the order of each of the following differential equations. 
& 
@) =8 () (B) +2 43-0 
Ce ay dy 
&) t=O (d) 253 tere 7 


By ‘an analytical solution’ we 
mean a solution which can be 
written in terms of a formula. 


If we allow p to be the zero 
function then Equation (4) 
includes the first-order linear 
differential equations 
considered in Unit 2. In this 
unit we wish to confine our 
attention to second-order 
equations and so we impose 
the condition that p is not the 
zero function. 


Obviously, letters other than y 
and x may also be used for 
the variables. 


6 MST204 6 


(ii) | Which of the equations in (i) are linear? 


(iii) Which of the following second-order linear differential equations are constant- 
coefficient equations? 


@ ry 
(a) xiatxty=0 © +4y=0 

dy ody : dy dy 
(0) 22 +3 + tins (d) 25a + 08x + Say =e 


(iv) Which of the equations in (iii) are homogeneous? 


[Solutions on p.49] 


Comment on theorems and proofs 

Certain important results about differential equations are stated in this unit as 
theorems. Since the objective of this unit is to teach you how to solve certain 
differential equations, proofs of theorems are not given. Theorems are stated where 
they give important information about the solutions of differential equations, and I 
will try to illustrate their importance. 


Many theorems require the condition that all the functions involved are twice 
differentiable throughout a suitable interval, and that this second derivative is 
continuous. I shall assume throughout the unit that all functions mentioned satisfy 
these conditions, without stating this explicitly in the theorems. 


Providing that you can solve the differential equations, understanding the 
theorems is only a secondary objective of this unit. No knowledge of the proofs is 
required. All the proofs omitted may be found in text books such as An 
Introduction to Linear Analysis by Kreider, Kuler, Ostberg and Perkins (Addison- 
Wesley 1966), and many also in M203: Introduction to Pure Mathematics. 


Study guide 


The first four sections of this unit are concerned with linear constant-coefficient 
equations. These are equations of the form 

dy dy 
ge tba ty =I) 
where a, b, and c are constants (and a ¥ 0). In Section 1 we examine the 
homogeneous case, where f is the zero function. The inhomogeneous case, where f is 
some non-zero function, is considered in Section 2. 


a 


Section 3 is concerned with finding particular solutions that fit given conditions 
d 
(such as y = 2 and ° = | when x = 0). There is an audio-tape associated with 


this section. Section 4 is the television section, and examines the graphs of some of 
the functions that arise as solutions of linear constant-coefficient equations. Section 
5 concerns a numerical method for second-order differential equations. Section 6 
contains additional exercises. 


The television section makes some reference to the material in Sections 1 to 3 and 
so it would help to have studied these sections before watching the programme. 
However, if you have not finished the first three sections by the time the television 
Programme is due to be broadcast, make sure you have read the preliminary notes 
at the beginning of Section 4 before watching the programme. The audio-tape 
cannot be studied until you have completed Sections | and 2. Although there are 
five sections of expository text you should not find the unit excessively long since 
Sections 3, 4 and 5 are all fairly short. 
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1 Homogeneous equations 


1.0 Introduction 
In this section we are concerned with the solution of equations of the form 
2, di ate general 
d jomogeneous 
a +b +ey=0 (1) constant-coefficient 
Oo dx second-order 
where a, b and c are constants, and a is not zero. I shall first describe (in differential equation. 


Subsection 1.1) the method whereby the general solution of an equation of this 
form can be found. Afterwards, in Subsection 1.2, I shall explain why the method 
works. 


1.1 The method of solution 


The auxiliary equation 
To indicate the idea of the method we consider the linear first-order differential 
equation 


dy 
a =O. 


One solution of this equation is y = e**. This suggests that we might try the same 
exponential form for the solution of the second-order Equation (1). 


To try y = e** in Equation (1) we first calculate 


dy, 
ae 
and 
2, 
oy = Be, 
so that 
dy dy a2 okx 5 ax 
apa tba t+ va aie + bie** + ce 


= (ai? + bA + cle*. 
Thus, y = e** is a solution of Equation (1) provided that 
ai? + bh +c=0. 


This equation plays such an important role in solving the differential equation that 
we give it a special name: 


Definition 1 
The quadratic equation 
ai? +bA+c=0 
is called the auxiliary equation of the Differential equation (1). 


Example 1 
The auxiliary equation of the differential equation 
dy dy 
322g + 47 =0 
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Solving the differential equation 

You should recall that the solution of a quadratic equation with real coefficients 
(such as the auxiliary equation) may be divided into three cases: 

(i) The equation has two distinct, real, solutions. 

(ii) The equation has one real solution (two equal roots). 

(iii) The equation has two conjugate complex solutions. 

These three cases for the solution of the auxiliary equation correspond to three 
cases for the solution of the Differential Equation (1). In the procedure below I 
describe how to write down the solution of Equation (1) in each case. Later in this 


section (Subsection 1.2) I shall explain why the solutions take the form they do, 
For now, though, I am only concerned with the technique. 


Procedure 1.1 
To solve the homogeneous differential equation 

oe a 5: +cy=0 

dx? Pax * = 

where a, b and c are constants, and a # 0. 
1. Write down the auxiliary equation 

a? +bh+c=0 
and solve it. 
2. The solution of the differential equation takes a different form in 
each of three cases, depending on the type of solution of the 
auxiliary equation. 


(i) If the auxiliary equation has two, distinct, real solutions 
4=A, and A= A,, then the general solution of the 
differential equation is 

y = Ae™* + Be. 

(ii) If the auxiliary equation has equal roots, so that it has just 
one solution 2 = a, then the general solution of the 
differential equation is 

y=(A+ Bxje* 

(iii) If the auxiliary equation has complex roots, A = « + if and 
4 =a — if, then the general solution of the differential 
equation is 

y = &*(A cos Bx + B sin fx). 
In each case A and B are arbitrary constants. 


The remainder of this subsection consists of worked examples illustrating the 
procedure above. You may wish to read through these, or you may prefer to try 
the examples yourself first, and refer to their solutions only if you get stuck. 
Example 2 

Solve the differential equation 


Method 
The auxiliary equation is 


#-314+2=0. 


Its solutions are A = 1, A = 2. This is Case (i) of Procedure 1.1 and so the general 
solution of the differential equation is 


y = Ae + Be* 
where A and B are arbitrary constants. 
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Example 3 

Solve the differential equation 
a 

Pe +4y=0. 


Method 
The auxiliary equation is 
#+4=0. 


The solutions of this are A = 2i, A = —2i. This is Case (iii) of the Procedure 1.1, 
with « = 0, 8 = 2. Hence the general solution of the differential equation is 


y =A cos 2x + B sin 2x 
where A and B are arbitrary constants. 


Example 4 
Solve the differential equation 
Foil ui du 


atts +u=0. 


Method 
The letters used for the variables are different here, but this makes no difference to 
the method. The auxiliary equation is 


47 +444+1=0. 


The left-hand side is the perfect square (24 + 1)? and so the only solution is 
A= —4. This is Case (ii) of Procedure 1.1, so the general solution of the 
differential equation is 


u=(A + Bre", 
where A and B are arbitrary constants. 


Example 5 
Solve the differential equation 
al z 
ae 3% at 2z=0. 
Method 
The auxiliary equation is 
247 -344+2=0. 
Using the formula for solving the quadratic equation gives 
3+ /9= 16 Se 16 3, v7 
ae aca 


This is Case (iii) of Procedure 1.1, with a = Dee, the general solution of 


the differential equation is 
3 JI tt 
z=exp (4) (400s (4) + Bsin ” 
where A and B are arbitrary constants. 


Example 6 (Where 2 = 0 is a solution of the auxiliary equation.) 
Solve 


@y od 
atig= 


You may wish to check the 
solutions obtained in these 
examples by differentiation 
and substitution. 
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Method 
The solutions of the auxiliary equation 2? + 34 = 0 are 4 = 0, A = —3. This is 
Case (i) of Procedure 1.1. The only difficulty lies in realizing that e* = e® = 1, so 
that the general solution is 

y=A+Be*, 
where A and B are arbitrary constants. 


Exercise 1 (Practice in solving homogeneous differential equations.) 
Solve each of the differential equations below. 


P ay 
a) 16% + 5y=0 w) 42% by =0 
(i) Gar +90=0 7) aratéee 
dz ’ 
(ii) #420 wi 223 43% a0 


[Solution on p.49)] 


Sometimes we need to solve equations which have coefficients expressed in terms 
of unspecified constants rather than numbers. This does not affect the method of 
solution, 


Example 7 

Find the general solution of the differential equation 
dx dx 
e + 40— = +o*x=0 


where @ is a constant. 
Method 
The auxiliary equation is 
2? + 404 + wo? = 0. 
The values of A aginst this are 


an ee 2 (-2 Vo. 


We have distinct, real, roots (Case (i) of Procedure 1.1). The general solution of 
the differential equation is therefore 


x = Ad-2+V3)en 4 Bel 2- JB, 


where A and B are arbitrary constants. 


Exercise 2 

Find is general solution of each of the differential equations below. 
dy 

0) i an 20 +o7y=0 


(ii) os T+ 60% + Sux =0 


where « is a constant. 
[Solution on p. 49] 


1.2 Why the method works 


We now look at why the expressions given in Procedure 1.1 for the general 
solution of a homogeneous equation (Equation (1)) are correct. Before looking at 
the specific expressions given for the solutions, I will discuss some ideas relevant to 
the solution of linear differential equations in general. 
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Definition 2 
If u, and w, are functions, then a linear combination of u, and u, 
is any function u of the form 


u(x) = c,u,(x) + c,u,(x) 


where c, and c, are constants. 


Example 8 
Suppose u,(x) = x, u(x) = sin x. The functions 
u(x) = 3x? + 2 sin x, 
u(x) = 2x? — 4 sin x, 
are linear combinations of u, and u,. On the other hand the functions 
u(x) = x? sin x, 
u(x) = x? + sin 2x, 
u(x) = 2sinx + x*, 


are not linear combinations of u, and u,. 


Theorem 1 


Suppose that the functions u, and u, are solutions of a This theorem is fairly easy to 
homogeneous linear differential equation. Then any linear oe 

combination of u, and u, is also a solution of this differential 

equation. 


The next example illustrates how this theorem can help us validate the method in 


Procedure 1.1. 
Example 9 
Consider the differential equation 
u"(x) — 3u'(x) + 2u(x) = 0. (2) Bomanber that u(x) means 
Each of the functions u,(x) = e* and u,(x) = e** is a solution of this differential 7 is 


equation. (This follows from the fact that 4 = 1 and 4 = 2 are solutions of the 
auxiliary equation 4? — 32 + 2 = 0.) We now know, from Theorem 1, that any 
linear combination of u, and u, is also a solution of Equation (2). That is, any 
function 

u(x) = Au,(x) + Buz(x) = Ae* + Be**, 
where A and B are constants, also satisfies Equation (2). For example, one 
solution is 2e* + 3e”*. 


In each of the three cases described in Procedure 1.1, the formula given for the 
solution of Equation (1) is a linear combination of two functions. 

Case (i): Ae*™ + Be. 

Case (ii): Ae* + Bxe*. 

Case (iii): Ae“cos Bx + Be™* sin Bx. 
Theorem 1 shows that so long as we know that the individual functions (e**, e*** 
in Case (i); e* and xe* in Case (ii); e“sin Bx and e*cos Bx in Case (iii)) are 
solutions of Equation (1), then these linear combinations must also be solutions. 
To check that the individual functions are in fact solutions is not too complicated. 
Let us look at this now. 


Checking the solutions 
The general homogeneous constant-coefficient linear second-order differential 
equation (1) may alternatively be written as 


au"(x) + bu'(x) + cu(x) =0 (3) 
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(a notation which emphasizes that the solutions u(x) are functions). Now we saw 
in Subsection 1.1 that u(x) = e** is a solution of Equation (3) so long as A satisfies 
the auxiliary equation 

a? +bh+0=0. (4) 
Let us now check that the functions specified in Procedure 1.1 are solutions of 
Equation (3) in the various cases. 
Case (i) 
In this case A = 4, and 4 = A, are real distinct roots of Equation (4). Thus the 
functions u,(x) = e** and u,(x) = e”* are already known to be solutions of 
Equation (3). 
Case (ii) 
In this case A = « is the only solution of Equation (4). Thus u,(x) = e* is known 
to satisfy Equation (3). The expression given in Procedure 1.1 for the general 
solution of Equation (3) also contains the function u(x) = xe*. We must check 
that this is also a solution. We have 

u,(x) = xe* 

u(x) = &* + axe™ = (1 + ax)e* 

u3(x) = we™ + (1 + ax)ae™ = (20 + a? xje™. 
Hence 

au3(x) + bu (x) + cu,(x) = e*(a(2x + a?x) + b(1 + ax) + ex) 

= €*((2ax + b) + (aa? + ba + c)x). (5) 

Now a is a solution of Equation (4); so ax? + ba +c =0. 
Also, we know that the auxiliary equation has coincident roots, so that 


b? — 4ac = 0. Hence its solution is a = — 2 and so 20x +b =0. 

Thus the right-hand side of Equation (5) above is indeed ‘zero, and so u, is a 
solution of Equation (3). 

Case (iii) 

We know that A=a+if and 
so the complex functions 


(x) = eet 


4=a—if are solutions of Equation (4) and 


and 

u,(x) =e 
are solutions of Equation (3). In this course we will normally restrict our attention 
to real solutions of differential equations, and so the complex solutions u, and u, 


are not suitable for our purposes. However, using Euler’s formula (see Unit 5, 
Section 4) we can write: 


1 1 
e&*cos Bx = zm) + piatx) 


and 


t 1 1 
e*sin Bx = 3) - 3 M2). 
Hence, by Theorem 1, the (real) functions e**cos 8x and e**sin Bx are also solutions 
of Equation (3). These are the functions appearing in the formula in Procedure 1.1 


The general solution 

So far we have seen why the functions given by Procedure 1.1 are solutions of 
Equation (3). To show that the procedure gives the general solution of Equation (3) 
we must show that there are no other (real) solutions. This is the consequence 

of another general theorem. 


To make this argument 
precise we should define a 


when z is a complex variable. 
If w(t) = u(t) + iv(t), where 
u(t) and v(t) are real functions 
of t, then we define 
dt) _ du) dt) 

dt de at 
It is then possible to show 
that many of the usual rules 
of differentiation still apply. In 
particular 

ip 

= = ie”. 
Furthermore, having defined 
the derivative of a complex 
variable, it 
meaningful to talk about 
complex solutions of 
differential equations. 


We shall normally invoke 
Theorem 1 in the context of 
real solutions, but the theorem 
is true for complex solutions 
as well. 
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The basic idea of the theorem is that once we have found two ‘genuinely different’ 
solutions of a linear second-order homogeneous differential equation, then the 
general solution is given by an arbitrary linear combination of them. To state the 
Theorem we need the following definition of what ‘genuinely different’ means. 


Definition 3 
Two functions u, and u, having the same domain are linearly uo eee Ba sa constant 
ae te = multi es 
independent if neither function is a constant multiple of the other. pute ate the aa 
domain and if there is a 
constant A such that 
Example 10 u(x) = Au, (x) 
(i) If uy(x)=e* and u,(x)=e™* then u, and w, are linearly independent. for all x in their common 
For there is no constant A such that e** = Ae* (or e* = Ae?*) for all real domain 
numbers x. 


(ii) If u,(x)=sinx and u,(x)=2sinx then u, and u, are not linearly 
independent since u(x) = 2u,(x). 


Theorem 2 
Suppose that D is an interval of real numbers and that the A definition of interval is given 
functions u, and u, both have domain D and are linearly in the study comments for 


Unit 6 in Section 5 of the 


independent solutions of the differential equation Preparatory Booklet. 


P(x)ul"(x) + g(x)u'(x) + r(x)u(x)=0 (for all x in D) (6) 
Suppose also that for every x in D, p(x) is not equal to zero. 


Then the general solution of Equation (6) is Au, + Bu, where A 
and B are arbitrary constants. 


I have stated the theorem in a form applicable to general second-order 
homogeneous linear differential equations. Notice that a second-order equation 
with constant coefficients inevitably satisfies the condition p(x) # 0 for all x. 
(If p(x) = a = 0, Equation (6) is not a second-order equation!) 


My main aim in this section has been to explain why the method given in 
Subsection 1.1 for constant coefficient equations works. However, the theorems 
are informative for non-constant-coefficient equations also. 

Example 11 

The general solution of 


x7u"(x) + xu'(x) — u(x) =0 (x > 0) (7) 
is 

Wx) = Ax +3 (x > 0). 
Explanation 


Let u,(x) = x and u,(x) = 1 be functions which have the same domain (x > 0). 


We can show that these functions are solutions of Equation (7) by direct 
substitution. (I leave you to do this.) 


Now u, and u, are linearly independent functions. Also, Equation (7) is of the 
form described in Theorem 2 (for in this case D is the interval (x > 0) and so 
p(x) = x? is non-zero for all x in D. 


Hence, by Theorem 2, the general solution of Equation (7) is 
ux) =4 + Bx (x > 0), 


where A and B are arbitrary constants. 
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Exercise 3 
Find the general solution of the differential equation 
ay dy 
2a xt y= 
aie Rt 0 (x >0) 


Hint: first look for solutions of the form y = x". 
[Solution on p. 49] 


Summary of Section 1 
The general solution of 
ay 


dy 
a toetay=0 


depends on the nature of the solution of its auxiliary equation 
a? + bh+c=0. 
The following table summarizes the method: 


General solution of differential 
equation. (In each case A and B 


Type of solution of auxiliary equation are arbitrary constants.) 


Two real roots: 4 = 4, and A =A, 
Equal roots (one real solution): 4 = « 
Complex roots: 2 = « + if and 4 =a — ip 


y = Ae 4 Bex 
y=(A + Bxje* 
y = e&*(Acos Bx + Bsin Bx) 


This method may be justified by: 


1. verifying that the various functions (e*'*, e***; e*, xe**; e™* cos Bx, e* sin Bx) are 
indeed solutions in each case; 


2. referring to Theorems 1 and 2 stated in the text. 


The concepts of linear combination (Definition 2) and linear independence 
(Definition 3) for two functions were defined in the text. 


End of section exercises 


Exercise 4 
(i) What is the auxiliary equation of the differential equation 
dy dy 
Lm cs bar =0. 


(ii) Solve this auxiliary equation. 
(iii) Write down the general solution of the differential equation. 
[Solution on p. 49] 


Exercise 5 

Find the general solution of each of the differential equations below. 
, @y ody a GZ dy 
(i) ge t?A tno a) t= 4 
ay ey 70 

(i) Fa - 16 =0 (iv) Gz =O. 


[Solution on p. 50} 
Exercise 6 
For which values of the constant k does the differential equation below have a general 
solution that involves sines and cosines? 
dy 
grt ke + 4y=0. 


[Solution on p. 50} 
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Exercise 7 
Given that y=x and y= x? are solutions of 
ofa 4 2y=0 (for all real x) 
is it legitimate to use Theorem 2 to write down the general solution of this differential 
equation? 
[Solution on p. 50} 


2 Inhomogeneous equations 


2.0 Introduction 
In this section we study the solution of equations of the form 
au"(x) + bu'(x) + cu(x) = f(x) (1) 


where a, b and c are constants. These differ from the homogeneous equations 
studied in Section | in that a function f(x) has been introduced on the right-hand 
side. 


2.1 The method of solution: ‘particular solution 

plus complementary function’ 
In essence, the method of solution is similar to the method used in Unit 1 for 
solving inhomogeneous recurrence relations. The method depends on our being 
able to find a particular solution of Equation (1). Adding any such particular 
solution to the general solution of the homogeneous equation 

au"(x) + bu'(x) + cu(x) =0 

(obtained by replacing f(x) in Equation (1) by zero) gives the general solution of 
Equation (1). 
Before examining the theoretical justification for this method I will work through 
an example that illustrates how it works in practice. 


Example 1 
Find the general solution of the equation 

u(x) — 3u'(x) + 2u(x) = 6. (2) 
Method 


Step 1: We first solve the homogeneous equation obtained by replacing the right- 
hand side of Equation (2) by zero. That is, we find the general solution of 


u"(x) — 3u'(x) + 2u(x) = 0. 


Apart from a change of notation this equation is the same as the equation in 
Example 2 of Section 1. There we obtained the solution u = u,, where 


u,(x) = Ae* + Be**, 
(A and B are arbitrary constants.) 


Step 2: Next we find any one solution of the Inhomogeneous equation (2). This is 
not always easy, but in this case it is quite easily seen that the constant function 


u(x) = 3 


is a solution. (For if u,(x) = 3, then u,(x) = 0 and u;(x) = 0. Thus 
u(x) — 3uj(x) + 2u,(x) = 6, as required.) 


Step 3: Finally we add the particular solution found in Step 2 to the general 
solution of the homogeneous equation found in Step 1. The result u,(x) + u,(x) is 
the general solution of the original Innomogeneous equation (2). Thus 


u(x) = Ae + Be* +3 
is the general solution of Equation (2). 
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The theoretical justification of this method relies on a theorem (Theorem 1 below). 
Although we are concerned with constant coefficient second-order equations in 
this unit, the method works much more generally, and so I will state the theorem 
for general second-order linear equations. To do so I will need the following 
definition. 


f 


Definition 1 
Let 
P(x)u"(x) + q(x)u'(x) + r(x) u(x) = fix) (3) 


be any inhomogeneous linear second-order differential equation. 
Then its associated homogeneous equation is 
p(x)ul"(x) + q(x)u'(x) + r(x)u(x) = 0. 


Theorem 1 

Suppose that u, is any solution of Equation (3) and that u, is the 
general solution of the associated homogeneous equation. Then 
u, + u, is the general solution of Equation (3). | 


We often need to refer to ‘the general solution of the associated homogeneous 
equation’ and so it is convenient to give it a special name: 


Definition 2 


The general solution of the associated homogeneous equation is 
[ known as the complementary function. 


With this definition we can make a succinct statement of the method for solving 
inhomogeneous equations of the form of Equation (1): 


Procedure 2.1 a | 


To solve the inhomogeneous differential equation 

au"(x) + bu'(x) + cu(x) = f(x). 
1. Find the complementary function using Procedure 1.1. 
2. Find a particular solution. 


3. General solution 
= Particular solution + Complementary function. 


Example 2 
Find the general solution of the equation 

u"(x) + 9u(x) = 9x. (4) 
Method 


The associated homogeneous equation is 
u"(x) + 9u(x) =0 
which has the general solution u = u,, where 
u,(x) = A cos 3x + Bsin 3x. 
This is the complementary function of Equation (4). 
A particular solution of Equation (4) is the function 
u(x) = x. 


(For, if u,(x) = x, then u(x) = 1, so up(x) = 0. Thus u;(x) + 9u,(x) = 9x, as 
required.) 
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The general solution of Equation (4) is, therefore, 
u(x) = x + Acos 3x + Bsin 3x, 
where A and B are arbitrary constants. 


The finding of the complementary function is a question of applying the method 
from Section 1. However the finding of a particular solution is another matter. In 
the two examples above, I have carefully chosen equations for which a particular 
solution can be found ‘by inspection’. This sort of approach is not going to work 
in general, though. The next subsection looks at the finding of particular solutions 
in a more systematic way. 


Exercise 1 

Consider the differential equation 
@ 

re +4y =8. 


(i) Write down the associated homogeneous equation and the complementary function. 
Find a particular solution of the form y = constant. 

(iii) Write down the general solution of the differential equation. 

(Solution on p. 50) 


There is one point that sometimes causes confusion. An inhomogeneous equation 
has more than one particular solution (of course!), Does it matter which one we 

choose to add to the complementary function? It does not; the general solutions 
obtained may look different, but are actually equivalent. 


2.2 Finding a particular solution 

I shall not describe a general procedure for finding a particular solution of 
Equation (1). The form of the particular solution depends on the form of the 
function f(x) on the right-hand side of the equation. I shall describe the sort of 
functions that work for certain specific forms of f(x). On the whole, the approach 
is ‘try something similar to f(x)’. 1 will start with a fairly simple case. 


Sfislinear (f(x) =kx +1 

The first example I will consider has f(x) = 4x + 2. The basis of the method is to 
try a particular solution which is similar to f. In this case f is linear so I will try an 
arbitrary linear function y = mx + n, and find (if possible) suitable values for the 
constants m and n for which this is a solution. 


Example 3 
Find a particular solution of the equation 
Py dy 
3 a Fae +2 
Solution 


We try a solution of the form 
y=mx+n 


where m and n are constants to be chosen so that the differential equation is 
satisfied. As a preliminary to substituting this function into the differential 


; dy dy 
equation we note that = mand a 0; hence 


Py dy 
Ba — gy + = 0— m+ (mx +n) 
= mx + (n — 2m). 


For the function mx + (n — 2m) to be equal to the function 4x + 2 on the right- 
hand side of the differential equation, we must have 


mx + (n — 2m) = 4x +2 (for all x). 
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This can only be true if 
m=4 

and 
n—2m=2. 


We can solve these simultaneous equations (by substituting m = 4 into 
n — 2m = 2) to give m= 4, n= 10. Thus 


y=4x+10 
is a particular solution of the given differential equation. 


It is always sensible to check calculations such as these by substituting back into 
the differential equation as follows. 


Check: If y = 4x + 10 then 


Py dy 
Sat 2g FY = 3x 0-2 x 4+ (4x + 10) 
=4x+2, 
as required. 
Exercise 2 


Find and check particular solutions of the form y = mx + n for each of the differential 
equations below 
a a ree Co ee 
(i) Fa ay ta e+ 3. Gi) Fat iy tye 
[Solution on p. 50) 

(Notice that in part (ii) of the above exercise you still need to look for a solution 
of the form y = mx + n. Although the function f(x) is just a multiple of x in that 
case, there is no solution of the form y = mx.) 


fan exponential (f(x) = ke**) 
In the next example, f(x) = 3e**. ‘Try something similar’ to f(x) again works. This 
time we try an exponential, y = me**, and determine the value of the constant m. 
Example 4 
Find a particular solution of 
dy 2x 
qa ty=se 
Solution 
We try a solution of the form 


y=me™, 


‘en as 2 59 2Y 2 set 2 
Then 7 = 2me » and re = 4me + 80 sa + y = 4me™* + me = Sme**. 


Thus y = me* is a solution of the differential equation if Sme?* = 3e?*, that is if 


m= 2 Hence 


y= je 


is a particular solution of the given differential equation. 


Exercise 3 
Find a particular solution of the equation 
dy 4dy = 
a =z FS Ey = Ze 


[Solution on p. 50} 


We try an exponential with 
the same exponent as f- 
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Sf sinusoidal (f(t) = k cos wt + sin wt) 

The most important case so far as applications are concerned arises when the 
function f is a sine or a cosine or a linear combination of them both. In the next 
example fis given by f(r) = 2sin3r. Although there is no cosine term we look for a 
particular solution of general sinusoidal type y = mcos 3t + nsin 31. 


Example 5 
Find a particular solution of 

dy d 

a + z + 2y =2sin3r. (5) 
Solution 
This time we try 

y = mcos 3t + nsin 3t. 
Then 

dy . 

— = 3ncos 3t — 3msin 3¢ 

dt 
and 

2. 

oe —9mcos 3t — 9nsin3t 
so 

2 

% +94 2y = (—7m + 3n)cos 3t + (—3m — 7n)sin 3t. 
We require this to equal 2sin 32, So we need 

—7Tm + 3n=0 

—3m — 7n =2. 


These equations give m = -3 and n= -s So a particular solution is 


3 en 

y= 55 083 - 99 5k. 
Before asking you to do any examples, I will show you an alternative approach for 
this example. This uses complex numbers, and is easier so long as you are 
confident with their use. The method is based on the idea that any sinuosoidal 
function can be represented by a phasor (see Unit 5, Subsection 4.3). We look for 
a particular solution of Equation (5) which is a sinusoidal function (as in the 
solution above). This time, however, we do the calculations using phasors. 
Remember that the phasor of the sinusoidal function 

y =mcos wt + nsin wt 


is the complex number z for which 


y = Re(ze) 
and is given by 
z=m-—in. 


With this in mind, we can look for a solution of Equation (5) of the form 
y = Re(ze™) (6) 


where z is a complex number to be chosen so that the differential equation is 
satisfied. We have 
dy 


ka! A ize") 
at Re (3ize™") 


See Unit 5, Subsection 4.3, 
Example 4. 
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and 
£Y Re (—9ze**). 


So (6) is a solution of Equation (5) so long as 
Re((—9 + 3i + 2)ze™*) = 2sin3¢ 


Now the phasor of 2sin 3t is —2i so if we write the right-hand side in phasor form 
we have 


Re((—9 + 3i + 2)ze**) = Re(—2ie**). 
This holds if 


(—7 + 3i)z = —2i 


With this value of z, Equation (6) gives a particular solution of the Differential 
Equation ©) Tus particular solution is a sinusoidal function with frequency 3 and 


phasor — z + x It can therefore be written 
3 7 
[pa —=sin 3t. 7 
y cos 3t 39 03 (7) 


This alternative approach may look a little complicated at first but in practice we 
can abbreviate the calculation considerably. You will find that you can get from 
Equation (5) to Equation (7) fairly easily by applying the following general 
procedure: 


a 
Procedure 2.2(a) 


To find a particular solution of 


ft +o + cy = kcos ot + Isin wt. 


1, Write the right-hand side in phasor form 
Re ((k — ile™). 

2. Try a solution in the form 
y = Re(ze). 

3. Substitute in the differential equation and simplify to obtain 
(—wa + iwb + c)z =k — il, 


and solve for z. 


4. Then the particular solution is the sinusoid with frequency 
and phasor z. That is, if z = X + i¥, the solution is 


y = Xcosmx — Ysinwx. 


1 


Although the exercises in this unit can be done using either of the methods I have 
described for Example 5, I recommend that you get used to the complex number 
method. It is always quicker, and sometimes much quicker (for example, for the 
sort of problems discussed in Section 4). 


Exercise 4 

Finda a ae solution of the equation 
dy 
awe +34 39 — 90531 


[Solution on p. 50} 
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Exercise 5 (Optional, additional practice.) 
Find a particular solution of the equation 
dy 
<2 tye in 2¢, 
ae Ly = 8cos 2t + 16sin 2¢. 


[Solution on p. 51] 
In this subsection we have looked at the problem of finding a particular solution 


of Equation (1) for certain specific forms of f(x). We have seen that the following 
procedure usually works 


Procedure 2.2 
To find a particular solution of 


ed ff ot + cy =fix): 
1. If f(x) = kx +1, try 
y=mx+n. 
2. If f(x) = ke*, try 
y=me™, 
3. If f(x) = kcos wx + Isin wx, try 
y =mcoswx + nsin wx, 
or replace the right-hand side by Re((k — il)e’**) and try 
y = Re(ze™) 
as described in Procedure 2.2(a). 


I pointedly say that this procedure ‘usually’ works, for reasons we shall now 
examine. 


2.3 Exceptional cases for particular solutions 


Occasionally, the functions suggested in Procedure 2.2 fail to provide a particular 
solution. Let us see an example of this difficulty, and how it may be resolved. 


Example 6 

Find a particular solution of the equation 
a 
a — 4y = 2c, 

Solution 


According to Procedure 2.2 the obvious guess is y = me**. If we try this, 
COE riper es Ce ernsnecy Oe 
then 7 = 2me and 773 = 4me’ *, giving 


ay 2x as 

me 4y = 4me** — 4me** = 0. 

So there is no value of m which makes y = me** a solution of the differential 

equation. What has gone wrong? The trouble is that the solution suggested in 

Procedure 2.2 contains the function e”* which is a solution of the associated 
a& 

homogeneous equation ze — 4y =0. It follows that me?* is part of the 

complementary function and so cannot be a solution of the inhomogeneous 

equation. 
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When Procedure 2.2 fails in this way we can usually overcome the difficulty by 
multiplying the function suggested in the procedure by x. In this case we try 


y = mxe**, 
Then 
2 = me** + 2mxe** = m(1 + 2x)e?*, 
and 
dy 2x 2x 2x 
aa = 2me + m(1 + 2x)e** = 4m(1 + x)e*, 
so that now 


12, 
a — 4y = 4m(1 + x)e** — 4mxe** = 4me?*, 
This can be made equal to 2e** by choosing m = 5 which gives the particular 
solution 

y= pxe™ 
In this example we see that the approach described in Procedure 2.2 does not 
work, but that by ‘multiplying by x’ we can still find a particular solution. 


Let us now look at two other examples where similar difficulties arise, and how 
they may be resolved. 


Example 7 
Find a particular solution of the equation 

dy 5 

Pr + 9y = sin 3x. (8) 
Solution 


Following Procedure 2.2 we would try y = mcos 3x + nsin 3x. However, since 

both cos 3x and sin 3x are solutions of the associated homogeneous equation 
2. 

2 + 9y = 0, this would just give zero when put into the left-hand side of 


Equation (8). To overcome this difficulty we again multiply the function suggested 
in Procedure 2.2 by x. That is, we try 


y= x(mcos 3x + nsin 3x). 
We can simplify the calculation by expressing this in the phasor form 
y = Re (xze™*) 


where z = m — in is to be chosen so that the differential equation is satisfied. We 
have 
; # = Re((1 + 3ix)ze**) 

and 

dy i Bix 

rr io Re ((6i — 9x)ze**) 
80 

d’y 

dx? 
The right-hand side of Equation (8) can be written Re(—ie**), If this is to be 
equal to Re(6ize**) we require 


+ 9y = Re (6ize**), 


6iz = —i 
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that is 

ee 

=a 

Substituting this into y = Re(xze**) gives the particular solution 

y= - cos 3x. 
Example 8 
Find a particular solution of the equation 

dy ody 

zat car =2x+2 (9) 
Solution 


This time we have a linear function on the right-hand side, so the first step is to 
try 


y=mxtn (10) 


d 2. 
as suggested in Procedure 22. Ify = mx +n then $= m and 4 = 0 giving 


Thus, for (10) to be a solution of Equation (9), we require 2m = 2x + 2 (for all x). 


Unfortunately this is impossible since 2m is a constant where as 2x + 2 varies with 
x. This shows that the differential equation does not have any particular solutions 
of the form y = mx + n. (As with the two previous examples, Procedure 2.2 breaks 
down here because part of the function y = mx + n is a solution of the associated 
homogeneous equation. The part of the function which cauges the trouble is y = n, 
ad dy 

for this is a solution of the associated homogeneous equation a + Ee =0and 
therefore does not survive substitution into the left-hand side of Equation (9).) To 
overcome this difficulty, I will again multiply the solution suggested in Procedure 
2.2 by x and try 


y = mx? + nx. 
2 
Then 2 = 2mx +n and 22 = 2m, s0 that now 
dx dx? 
dy ody 
gett ag amet: 
For y = mx? + nx to be a solution we require 
2m + 2(2mx +n) = 2x + 2, 
that is 
4mx + (2m + 2n) = 2x + 2. 


This can only be true for all values of x if m and n satisfy the simultaneous 
equations 

4m =2 

2m + 2n = 2. 
Hence m = 4 and n = 4. Thus the differential equation has the particular solution 


ya30? +9. 


In each of the examples in this subsection the method described in Procedure 2.2 
did not work because the solution we tried contained a term that was a solution of 
the associated homogeneous equation. In each case we were able to find a solution 
using the procedure on the next page. 
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Procedure 2.3 

If the approach to finding a particular solution described in 
Procedure 2.2 fails to produce one, then try x times the function 
specified there. 


Exercise 6 


Find a particular solution of each of the differential equations below. 
dy 

ey re z 

i) ae + 4y = 2008 2. 


[Solution on p. 51] 
That is very nearly the whole story. The next exercise completes it. Whether or not 
you try the exercise yourself, do look at its solution. 
Exercise 7 (Challenge) 
Find a particular solution of the equation 
@y dy 
me a +y=e. 
[Solution on p. 51] 
2.4 The superposition principle 
There is one final result that is relevant to finding particular solutions. I shall state 


the theorem using a compact notation for the left-hand side of a general second- 
order linear differential equation. I will write 


a d 
LO) = pla) + ale) + roxy 


Theorem 2 
Suppose that 
L(y) = fix) 


and 


L(y) = g(x) 
are any two linear inhomogeneous differential equations, with the 
same associated homogeneous equation L(y) = 0. Suppose also 
that y = u(x) is a solution of the first equation and y = v(x) is a 
solution of the second equation. 
Then, if a and b are constants, y = au(x) + bv(x) is a solution of 
the differential equation 


L(y) = af(x) + bg(x). 


This result is known as the principle of superposition. It states that if the right-hand 
side of our differential equation is a linear combination of functions for which we 
know particular solutions, then the same linear combination of the particular 
solutions is what we want. 


Example 9 
Find a particular solution of the equation 
2 


dy : 
ae + 9Y = 4sin3x + 13. 


MST204 6.2 25 


Solution 
A particular solution of the equation 


. 
oY + 9y = sin3x 


isy=— J xcos3x (see Example 7), A particular solution of” 
2 


ay 
pet =i 


isy= i (Check it!) Then Theorem 2 tells us that 


aod 
orang 9 


is a solution of the given differential equation. 


When the right-hand side of an equation is the sum of two parts we can either 
solve for the two parts separately, as in Example 9, or together, as illustrated by 
the following example. 


Example 10 
To find a particular solution of the equation 
c art 3 4 we de™ — 78in 24, 
dt dt 
we would try 


y = le™ + moos 21 + nsin2r 
(or y = le™ + Re (ze), 


Exercise 8 

Find particular solutions of each of the differential equations below. 
1) oe oe 

(i) a -25 +y=4e—3e* — (Use the result of Exercise 7). 


(i) a 5 0% 5 100% = 9co0s4r + 100 


[Solution on p. 51} 


2.5 General solutions 


In Subsection 2.1 we discussed how the general solution of an inhomogeneous 
equation of the form given in Equation (1) on page 15 can be found by adding a 
particular solution to the complementary function. Now that we have developed a 
method for finding particular solutions, you should be able to find general 
solutions for suitable functions f(x). 


It is sensible when finding a general solution to find the complementary function 
first. This will provide advance warning of the sort of ‘exceptional case’ discussed 
in Subsection 2.3, 


Exercise 9 (Consolidation of all of Section 2) 
Find the general solution of each of the differential equations below. 


(i) oO 40 = 2t «iy ys aye 2e* — Se** 
iy oY PLY hey tv) 4 20 4 20%y = Ssin 2a + wh 


dx 
ioe @and ¢ are constants) 


[Solution on p. 52} 
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Summary of Section 2 

1, The general solution of an inhomogeneous linear differential equation is found 
by adding any one particular solution to the complementary function 
(Theorem 1). The complementary function is the general solution of the 
associated homogeneous equation. 

2. To find a particular solution of an inhomogeneous constant coefficient 
equation, we try a particular expression for y, and see whether it works. The 
form to try for y depends on the form of the function on the right-hand side of 
the equation. The table below gives expressions that usually work. 


kcos@x + Isinwx| {or 


Re (ze’”*) 


{ere + Pe 


3. In some cases there is no particular solution of the form suggested in the table. 
This occurs when the function to be tried contains a term that is a solution of 
the associated homogeneous equation. 


If the expression given in the table for y does not work, then try x times the 
given expression. 

4. If the right-hand side of the differential equation is a linear combination of 
functions given in the table above, then a particular solution may be found by 
trying the same linear combination of the expressions given there for y 
(Theorem 2). 


End of section exercises 


Exercise 10 


What form of expression for y should you use, in order to find a particular solution of each 
of the differential equations below? 


(i) £y ay mes, (iii) fy ay ae™, 

Gi) © — 4y = sind, Gv) 2 4 4y = sind + 3x 

[Solution on p. 53} 

Exercise 11 

Find the general solution of each of the differential equations below. 
dy ody 


ji) get tat 
ie : 
(i) Sy 4y =sin3x. 


[Solution on p. 53] 
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3 Initial and boundary conditions (Tape Section) 


In this section we will look at solutions of second-order differential equations 
satisfying additional conditions, In Unit 2 you saw that the imposition of an extra 
condition (such as y = 2 when x = 0) on the solutions of a first-order differential 
equation usually picks out one particular solution by sétting a value on the one 
arbitrary constant in the general solution. 


Now the general solution of second-order differential equation has two arbitrary 
constants. So to pick out a particular solution from the general solution we might 
expect to require two additional conditions. This is indeed the case. 

In the first part of the tape I will work through two examples which illustrate how 
we can find solutions satisfying given initial conditions. These examples will also 
provide you with an opportunity to review the techniques in Sections 1 and 2. In 
the second part of the tape I will discuss the uniqueness of solutions satisfying 
given conditions. 


Start the tape when you are ready. 
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( Example | 


Find the particular solution of 


2 
d 
¥ = grey 


+ 2y = 
dx* ax 2 


Which satishes y = 


Solution 


ausciliary 
equation 


general 
Solution 


Particular 
Solution 


(2 ] Example 2 


Find the particular solution of the inhemegeneous equation 


2, 
sone ay 
at* 


for which y =° 
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3 ] Inhomogeneous equations 


Find the complementary function. 
Find any particular solution. 
Add these to get the general solution. 


Then apply the conditions to obtain the required 
particular solution. 


Note: the particular Solutions obtained in steps 2 and 4+ 
are no€ necessarily the Same. 


The general selution of 
dt> 


Avaaliary 

equation is 

At +d associated 

ieee homogeneous 
equation 


ay Gus & 
A. ee = (© +e 


So we require L=%o, m= Yq and n=0 giving 


= e particuler 
Y= Hoe + %e Solution 


general 3 iat 
Solution YU=Acos B+Bsin3t rise tqe | 


| 
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rs ] Applying the conditions 


general solution Y= Acos 3t+Bsin at+ Noe + at 


O=A cosO+Bsino+ Ke + 4x0 


=A+%X% 


So the particular solution satisfying the conditions is 


7 
pe agi} + 
27o SIN Bt -75 cos BE tee + He 


nn 


from ¥ 
Complementary original 
function particular 
Solution 
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© crecking lution. 


PROBLEM: seve d'Y yqy=e +t; Y= 3 =) when t<o 
dt* 


7! : oy t 
SOLUTION : ¥ =2o Sin Bt -io coset +¥%e 


Checking the expression satisfies the equation. 


1. Sin 3t and cos 3t are solutions of 


\sS @ particular selution 
ete +r whe rht) = eket 
So the expression satisfies the differential equation 


check the solution satisfies the conditions 


pe sin 3t- 73 


a “ei FUSES ee 
x 3cos 3t+ 7G Sin Bt +iQe aa 


Now work through Exercise 1 and check your solutions to make sure that you have 
understood the ideas so far. 


Exercise 1 

In each case find the solution of the differential equation satisfying the given conditions. 
(i) u"(t)+4u(t)=0;  u(n/2)=0 and —sw'(n/2) = 1, 

., @y dy . dy 

@) 72-3 tas y= 4 and 7 = —1 when x = 0. 

[Solution on p. 53] 


Now try Exercise 2, which is an introduction to the next part of the tape. Its solution 
is covered in the tape commentary. 


Exercise 2 
Can you find a solution of 
u(x) + 4u(x) =0 


such that u(0) = 0, u(x/2) = 17 


After you have tried Exercise 2, restart the tape. 
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@ Answer to Exercise 2 


The general solution is u(¥) =Asin 2x +B cos 2 


© =A sino +B coso 1=A sin t+ & cos 


so 
& 


Hae ula) =k and w (a) =L 
Conditions Same valve of =< 


u(a)=k and ule’) = & 


or 


wla)=K and ui(b) = L 


initial Condition problems 


het “ ' 
P(x) w (x) +q4Gu (x) + r()uG = Ff) 


be any linear second-order differential equation, for 

which p(x) does not take the value 2ero for any valve 
of x. Then there is one and only one Solution of this 
differential equation satisfying the pair of initial Conditions: 


u(a) and u(a)=Lb 
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Exercise 3 
Identify each of the following as either an initial condition problem, or a boundary 
condition problem. Solve each problem. 

(i) u"(x) + 4u(x)=0; uO) =1, v0) = 0. 

(i) u"(x) + 4u(x)=0; — u(0)=0, —u(n/2) = 0. 

(ili) w"(x) + 4u(x)=0; — u(0)=0, —-w’(0) =0. 

(iv) u"(x) + 4u(xy=0; (=n) = 1, win/4) = 2. 


[Solution on p. 53] 


These examples bring up some points of interest. Boundary condition problems 
may be ‘well-behaved’, in that they produce a unique particular solution, as in (iv) 
above. Indeed, this is often the case. However, a boundary condition problem may 
have no solution, as we saw in the tape, or it may have infinitely many solutions, 
as in (ii) above. 


For a homogeneous equation, with initial conditions of the form u(a) = 0, u'(a)=0 
there is no need to do any detailed calculation, as the solution must be the zero 
function (c.f. example (iii) above. This short cut cannot be applied to boundary 
condition problems, as there is no guarantee that the zero function is the only 
solution (c.f. example (ii) above). Nor can it be applied to inhomogeneous 
equations, for then the zero function is not a solution at all. 


Exercise 4 

Solve each of the following problems, 

, ay dy 

() Fa-y= 0; y=Oand * =Oarx=0, 

ee R dy 

Gi) Gaya 8 y= Oand 7 =Oarx=0. 
2. 

(iii) oF 20 + ox <0; x=aand = oar=o, 
where w and a are given constants. (You can use Exercise 2(i) of Section | here, if 
you wish.) 


{Solution on p. 54) 


Exercise 5 (Challenge—optional) 


(i) There is an analogue of Theorem 1 for first-order equations. What would you expect 
it to say? 

(ii) Can you find any solution(s) of 

dy 

rs -y=0 
satisfying y = 0 when x = 0? 

(iii) How does the answer to (ii) relate to your ‘existence and uniqueness’ theorem in (i)? 

[Solution on p. 54) 


Summary of Section 3 


In this section we examined the problem of finding solutions which, in addition to 
satisfying a second-order differential equation, also satisfy two extra conditions. 


Such conditions are called initial conditions if they both occur at the same value of 
the independent variable. Otherwise they are called boundary conditions. 


For initial conditions we have an existence and uniqueness theorem which is 
stated in Frame 9. 


For boundary conditions there is no guarantee that a solution will exist. Even if a 
solution does exist, it is not necessarily unique. 
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End of section exercise 


Exercise 6 

Find the solutions (if any) of each of the following problems. 

(i) u(t) + 4u'(t) + Su(r) = 0; u0)=0, u'(0)=2. 

(ii) u(t) + 9u(t) = 0; u(0)=0, — u'(x/3) = 1. 

(iii) Hayao; y=0, ® omx=n 

State in each case whether it is an initial or a boundary condition problem. 
[Solution on p. 54) 


4 Oscillations and graphs (Television Section) 


4.0 Introduction 


In this section we study the graphs of certain functions that arise as solutions of 
second-order differential equations. We will concentrate on cases of importance in 
dynamics, particularly those where the solutions contain sinusoidal terms, and so 
are oscillatory. Many of the ideas discussed in the programme are also covered in 
the first three sections of the unit. Those ideas introduced in the programme, and 
not covered elsewhere in the unit, are summarized in Subsection 4.2. 


4.1 Preliminary work 
The programme divides into three parts. 
1, We study the differential equation 


2. 

es —ky=0. (1) 
If k is positive the solution of Equation (1) involves exponential functions, but if k 
is negative it is sinusoidal. The programme reviews some ideas associated with 
sinusoids that were introduced in Unit 5. It also examines the relation of the 
amplitude and phase of the solution to initial conditions specified on the solution 
and its derivative. 


2. We study the homogeneous equation 
dy 
ax? 

(The algebraic forms of the coefficients are chosen because they provide simplicity 

later.) We know that the solution of Equation (2) is oscillatory (in fact purely 

sinusoidal) if a = 0, because we then have Equation (1) with k negative. For 


0.< « <q the solution is an oscillation whose amplitude decreases with increasing 
x, but if a > @, the solution is not oscillatory. 


+ rat +o*y =0. Q) 


Either way, so long as a > 0, the solutions decrease to zero as x increases. 
3: We study an inhomogeneous equation: 
dy ody : 
ga t 4g + Sy = sin2x. (3) 


Any solution of Equation (3) is the sum of two parts. One part is sinusoidal with 
the same period as sin 2x, and is the same for all solutions of the equation. The 
other part is a solution of 


2 
1 4B syao 


This part is an oscillation whose amplitude decreases as x increases, and becomes 
very small as x becomes large. (This second part of the solution may therefore be 
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called transient.) For large x only the first part of the solution remains (and so this 
part is sometimes called the steady-state solution). 

To gain maximum benefit from the programme you will need to have studied the 
Preceding sections of the unit. If you have not, you will find many of the ideas 
from these sections reviewed in the programme, but you will probably find that 
the pace of the programme is too fast for you to follow the details. Either way, 
temember that the important things to concentrate on are the shapes of the graphs 
of the solutions of Equations (1) and (2), and the way a solution of Equation (3) 
can be divided into two parts. 

The programme makes some use of the ideas concerning sinusoidal oscillations 
and phasors that where introduced in Unit 5. Exercise | serves to remind you of 
these. Before viewing the programme, work through this, and, if you have studied 
Section 1 of the unit, Exercise 2 as well. 


Exercise 1 
Find the phasor of the sinusoidal function 


S(t) = 4cos St — 3sin St, 
and hence find the amplitude and phase of this oscillation. 
[Solution on p. 55) 


Exercise 2 

Solve the auxiliary equation of the homogeneous differential equation 
#Y oY sot 
Pr + uae +o7y=0. 

For what range of values of « and @ are the roots real? 


[Solution on p, 55) 

As soon as possible after viewing the programme, read the notes in Subsection 4.2 
and then work through Exercises 3, 4 and 5. The remaining exercises may be left 
till later if you wish. If you have not studied the earlier parts of the unit it is still 
advisable to study the notes and try Exercises 3, 4 and 5 after viewing the 
programme, although there will be points that you will need to return to after you 
have studied the first three sections of the unit. 


Now watch the television programme ‘Good vibrations’. 


4.2 Good vibrations 


This subsection summarizes the new material in the television programme that you 
will be expected to know, and also takes some of the ideas a little further. 


1. An equation that will be useful in Unit 7 is 

@ 
a +o%y=0. (4) 
The general solution of this equation (see Section 1) is 

y = Ccosmx + Dsinwx, 


where C and D are arbitrary constants. You saw in Unit 5, Subsection 4.3, that 
this can also be written as 

ye eee Os 6) 
Thus (5), with A and ¢ arbitrary constants, is an alternative way of writing the 
general solution of Equation (4). This solution is a sinusoidal oscillation (see 
Figure 1). The constant q in (5) is called the angular frequency of the sinusoid, and 
is determined by the differential equation itself. The amplitude A and phase ¢ of 


the oscillation depend on the initial conditions (the values of y and 2 at x =0). 
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“? 


Figure 1. The graph of the sinusoidal oscillation y = A cos (wt + ¢) 


Fitting different initial conditions will have the effect of expanding or contracting 
vertically the graph in Figure 1 (varying A), and shifting the curve sideways 
(varying ¢) as in Figure 2. 


Figure 2. The effect of varying the constants A or @ in Figure 1 


2. An equation that will be useful in Unit 8 is 


yd 
Sa tage + oty = 0. (6) 


The nature of the solutions of this equation depends on the signs of « and 
a? — w*. Of particular importance for dynamical applications is the case when « is 


positive. The fundamental nature of the solution in four possible cases is 
summarized below. 


negative 

See the solution to Exercise 12 
at the end of this section for an 
eras ition of the entries in this 
tabl 


The graphs of typical solutions in each of these four cases are shown in Figure 3. 
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”Y 
*Y 


(a) a>0, 2 -w?>0 (b) «<0,27-w?>0 


(c) 2 >0, 0? -w? <0 (d) «<0,2?—w? <0 


Figure 3. Graphs of various solutions of Equation (6) 


Let us concentrate for a moment on the case when « > 0, a? — w? > 0 

(Figure 3(a)). In this case the general solution of Equation (6) is a linear 
combination of two decreasing exponentials. By choosing different linear 
combinations we can find solutions fitting different initial conditions. These 
solutions have graphs which are variants of the graph in Figure 3(a). Figure 4 
shows the sort of graphs that we can obtain in this case, All the solutions have the 
property of tending to zero as x becomes large. 


ye ys 
x 
> 
x 
ye v4 
- 
x x 


Figure 4. Graphs of possible solutions of Equation (6) when 2> 0 and a? — w? > 0, fitting 
various initial conditions. 


In the cases where « <0, a? — «* > 0, we get solutions showing similar variations 
on the basic increasing exponential shape shown in Figure 3(b). In the oscillatory 
cases (Figure 3(c) and 3(d)), the general solution of Equation (6) can be written in 
the form y = Ae~* cos (vx + @) where v = \/@? — a? and A and @ are arbitrary 
constants. Changing the values of A and ¢ will alter the initial conditions 
satisfied by the solution, but the solution curve will have the same basic shape. 


38 MST204 6.4 


; Ae 


Figure 5. In oscillatory cases, changing the initial conditions changes the values of the 
arbitrary constants A and @. The period remains the same. 


3. We saw in the programme, how a solution of 


ay | dy 
det ae 
is the sum of two parts (see Figure 6). 


y 
— 
x 


(a) A solution of Equation (7) 


pa 


) Its transient part 


+ Sy = sin 2x (7) 


= 


af 


ay 


(c) Its steady state part 


Figure 6. The solution (a) is the sum of two parts, shown in (b) and (c) 
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One part is a decreasing oscillation (Figure 6(b)), the other is a pure sinusoid 
(Figure 6(c)). For sufficiently large x the solution becomes effectively equal to the 
sinusoidal part. Whatever initial conditions we apply to Equation (7), this 
sinusoidal part of the solution is exactly the same. (That is, its amplitude, phase 
and frequency do not depend on the initial conditions.) The initial conditions 
affect only the decreasing part of the solution. 


The decreasing part of the solution is called a transient. The sinusoid to which the 
solution settles when x is large is called the steady-state solution of Equation (7). 


These ideas apply more generally. The general solution of 


dy 4 dy F 

gat tq, + Oy = keos vx + Isin vx (8) 
is the sum of the complementary function—the solution of 

BY, 5 Bi a 

te + ar +o*y =0, (9) 


and any one particular solution of Equation (8) (see Section 2). We have seen that 
if « > O the solutions of Equation (9) tend to zero as x becomes large. In other 
words, if « > 0 the contribution from the complementary function is a transient, 
There is one sinusoidal particular solution of Equation (8), and this is the steady- 
state solution. Any other solution of Equation (8) is the sum of a transient (from 
the complementary function) and this steady-state solution. The steady-state 
solution does not depend on the initial conditions. Only the transient term is 
affected by the initial conditions. 


4. Since the steady-state solution of Equation (8) (with a > 0) is a sinusoidal 
function, it can be represented by a phasor. The easiest way to find the steady- 
State solution is to calculate this phasor directly. The procedure for doing this is 
discussed in Subsection 2.3, and also in Example 1 below. Example 1 also 
illustrates how the phasor can be used to find the amplitude and phase of the 
steady-state solution. The angular frequency of the steady-state solution of 
Equation (8) is v, the angular frequency of the sinusoid on the right-hand side of 
the equation. 


Example 1 
Find the amplitude, phase and angular frequency of the steady-state solution of 
dy dy 25 
—> +9 = =x. 
set ax + 100y 9cos 7. 


Use the results to write down the steady-state solution. 


Solution 
The steady-state solution is a sinusoid with the same angular frequency as the 
sinusoid on the right-hand side of the equation, that is 


25 
angular frequency = T 
3 ‘ 25 225 
We are therefore looking for a solution of the form C cos Fad + Dsin qo 
equivalently A cos x + @). Since we are interested in the amplitude and phase 


we want the solution in the second form here. To do this we shall find the phasor 
z of this sinusoid by trying y = Re(ze'?*"*), which leads (by Procedure 2.2(a)) to 


2 
(+ + of) + 100}: =9. 
Solving for z gives 


2 P 
= — 95 (1 +i). 
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The amplitude and phase of the oscillation are given by the amplitude and phase 
of this phasor, that is 


amplitude = |z| = ~— 


The gai solution is therefore 
2h os i — 3x 


—x — —-}. 
Exercises 3, 4 and 5 are to be worked directly after watching the programme. 


Exercise 3 
Describe (in one sentence) the graph of any solution of the differential equation 


[Solution on p. 55) 


Exercise 4 
(i) What conditions must « and q satisfy if the general solution of the equation 
a + aut +o7y=0 


is to be a decaying oscillation? 
(ii) Sketch the graphs of the solutions of each of the differential equations below, that 
ge the specified conditions. 


: dy 
@ie oo y 4and Oatx=0, 


= 


B42 ey 0; y=2and Y= 1 atx=0. 


(Do not solve the equations. Your sketches should show (J) behaviour near = 0, 
(ID) behaviour for large ¢, (II) whether or not the solution is oscillatory.) 


[Solution on p. 55} 


Exercise 5 
Describe the solutions of the equation 
@y dy 


az qa tm toy= COS Vx, 


where 0 <a <, In particular, comment on the long-term behaviour of the solution, and 
the nature of the part of the solution that is dependent on the initial conditions. 


[Solution on p. 55] 


The following exercises are to be worked after viewing the een: and after 
working the first four sections of the unit. 

Exercise 6 

Find the amplitude, phase and angular frequency of the particular solution of 


PS fue doe 0 


satisfying y = 4 and = 9 at x =0. 
[Solution on p. 55} 


Exercise 7 
Consider the differential equation 
u"(x) + w?u(x) = 
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with the boundary conditions u(a)=k and ule + ca = |. Do these two conditions 
guarantee a unique solution? (Consider the cases k = | and k # | separately.) 

[Solution on p. 55] 

Exercise 8 

Find the particular solution of 


Bx. wes 
we + ar +o°x=0 
ae ee, . dx 
satisfying the conditions x =@ and qe 7 Oat =O. 
(Use the general solution of this equation found in Example 7 of Section 1.) 
[Solution on p. 55] 
Exercise 9 
Find the particular solution of 
#x (a 
at por +o*x=0, 


where 0 < <1 and 0 <oy that satisfies x = a and “* = 0 at ¢=0. 


Sketch the graph of this solution. 
[Solution on p. 56) 
Exercise 10 
Find the amplitude, phase and angular frequency of the steady-state solution of 
@x dx 
bard + wy + 100x = 9cos St. 
[Solution on p. 56} 
Exercise 11 
Find the phasor of the steady-state solution of 
dx i, 
at + 2Bag, + wx = cos vt 
where f and o are positive constants, and f < 1. Find the amplitude of this solution. 
[Solution on p. 56] 


Exercise 12 (Optional) 
Explain why the graphs of the solutions of the differential equation 


LY on® 4 3 
de + ma +o7%y=0 
take the forms specified in the table on p. 36. 
[Solution on p. 56) 


4 


Compare this with Theorem 1 
of Section 3. 
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5 A numerical method 


5.0 Introduction 


All the second-order differential equations you have encountered in this unit so far 
have been amenable to analytical solution, in that we have been able to express 
the solution in terms of a formula. Unfortunately, not all second-order differential 
equations which occur in practice can be solved analytically. For example the 
differential equation 

dy 5 

7) ies 10sin y, 
which arises in the study of the large oscillations of a pendulum, has no solution 
in terms of the elementary functions such as trigonometric or exponential 
functions. 


In such a situation we may resort to using a numerical method. In Unit 2 we used 
Euler's method to solve first-order differential equations of the form 


oY missy) 


This involves using a recurrence relation to generate values y, that give 
approximations to y(x,), where the difference x, ,, — x, between successive values 
of x is a constant h (the step length). In Euler’s method we approximate the 


derivative of y at x = x, by Yee = ~+ to obtain the recurrence relation 


Yoo =, + hm(x,,y,)- 


In this section we will show how an extension of this method can be used to solve 
second-order differential equations. 


5.1 Euler’s method 


There are several methods of obtaining numerical approximations to the solution 
of a second-order differential equation. The method I will describe here is to 
replace the second-order equation by two first-order differential equations, and to 
apply Euler’s method to each of these. (Later in the course, we will see other, and 
better, methods than Euler's for first-order equations. Use of these, together with 
the basic idea of replacing a second-order equation by two first-order ones, will 
provide more satisfactory methods for second-order equations.) 


The key idea is to introduce another variable z to represent the first derivative of 
y. That is, we write 
dj 
z= >. 
If we now differentiate z, we obtain 
dx dx™ 
To see iow this idea is used we look at the following example. 


Example 1 

Replace the second-order differential equation 
@y ody 
ee ag HO 

by a pair of first-order equations of the form 
ay =z 


dx 
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and 

a = (some function of x, y and z not involving their derivatives) 
Solution 
We can rewrite the differential equation as 

dz dy 

Pre bess + 6y=0. 


Also, replacing a by z, we obtain 
dz 
oe 5z+ 6y =0. 


Then the pair of first-order equations (written in a form suitable for applying 
Euler's method) is 


dy 
ae =z 
dz 
57 = Sz — by. 
r= 5z — 6y. 
Similarly any linear differential equation of the form 


& 4 
Paar + alvdge + rixdy = fix) 


(for which p(x) is non-zero for all x) can be written in the form 


dy 
ax * () 
é = m(x, y, Z) 


where m(x, y, z) is some function of the three variables x, y and z not involving their 
derivatives. Many non-linear equations can also be treated in this way. 


Exercise 1 


Replace each of the following second-order differential equations by a pair of first-order 
equations of the form given in Equation (1). 

. ay dy 

(i) gat tay t AY = 08x 

2 

a + 32siny=0 


[Solution on p. 57) 


(ii) 


Having changed our second-order differential equation into a pair of first-order 
equations, we can set up the equations for Euler’s Method. 


Example 2 
Formulate Euler’s method for solving the equation 
@y dy 
Bet Y= 0 2) 
Solution 
From Example 1 we can replace the equation by the two first-order equations 
dy 
m7? (3) 
dz 


Rn. (4) 
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To apply Euler’s method to Equation (3), we approximate the derivative of y at 
x = x, by 


Yet — Vr 
c (5) 


Here y, is Euler's approximation to y(x,), the value of the true solution at x,. The 
values of x, are given by 


xX, =Xq + rh. 
Replacing x in Equation (3) by the Expression (5) gives the recurrence relation 


Yous =, + he, 

where z, is the approximation to z(x,). 

Treating Equation (4) in a similar way we obtain 
2,41 = 2, + h(Sz, — 6y,). 


Thus Euler's method applied to the Second-order equation (2) yields the pair of 
recurrence relations: 


Yr = Sp + he, 
2,41 = —Ghy, + (1 + Sh)z,. 


Exercise 2 
Bers the recurrence relations for Euler's method applied to the differential equations 
y 


(i) xt 
Se dS + cot 
= iy + cos x 
(ii) wes 
& = —22siny 


[Solutions on p. 57} 


We are now almost ready to carry out some computations. All that remains is to 
see how the initial conditions fit into this scheme. Although my main purpose in 
developing a numerical method is to be able to attack equations that cannot be 
solved analytically, I shall deliberately illustrate the method with an equation that 
we can solve analytically, We will then be able to see how accurate the method is. 


Example 3 


Use Euler’s Method to obtain the value of y when x = | using a step length 
h =0.1 for the differential equation 


with the initial conditions y = | and % = 1 when x =0. 


Method 
We saw in Example 2 that Euler’s method gives rise to the recurrence relations 


Yer = Vp + hz, 
2,41 = —6hy, + (1 + 5h)z,. 


The initial conditions can be written as yy = 1 and zp = 1. Thus with h = 0.1 we 
have 

Yori =), + 0.12, 

Z41 = —0.6y, + 1.52, 


(¥o = 1 and zy = 1). 
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We can now use the recurrence relations to calculate y, and z,, then y, and z,, 
and so on. 


The following table (in which the results have been calculated to 3 decimal places) 
shows the results of the calculation. 


0.6y, + 1.52, 


5.2 Accuracy of the method 
For the equation 
dy dy dy 
ee = y= and = 1 when x =0 
considered in the previous subsection we can compare our numerical results with 
the analytical solution, which is 


y = 20% — 


The following graph gives both the analytical and numerical results. 


rit 


rot 


HERE: E 


Figure 1 


These results are not particularly good for this value of h though it can be seen 
that the numerical results do follow the general shape of the solution. With smaller 
step sizes we would expect the results to be much closer to the true solution 
(provided we work to enough decimal places). However, the smaller the step-size 
the more work we have to do to evaluate y at x = 1. Thus in general we would 
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wish to use a computer to obtain our results. The following table gives the 
approximation to y(1) obtained for various step sizes h. 


—5.247 


The correct value for y at x = 1 is 2e? — e? = —5.307. 


Exercise 3 
(i) Continue the table below to find an approximation to y(3.2), where y is the solution 
of the differential equation 
@ 
x +y=0 


TF 
x ys i Ye +022, |—-02y, +2, 
0 0 0.2 10 
0.2 0.2 04 0.96 
04 04 0.592 0.88 
0.6 0.592 0.768 0.762 
08 0.768 0.920 0,608 
10 0.920, 1.042 0.424 
12 1,042 1.127 0.216 
14 1.127 1.170 —0.009 
16 1.170 1.168 —0.243 
18 1,168 L119 -0.477 
2.0 1.119 1.024 —0.701 
22 1.024 0.884 —0.906 
24 0.884 0.703 —1.083 
26 0.703 0.486 
28 0.486 
3.0 
3.2 


(ii) What is the exact solution of this equation satisfying the given initial conditions? 


(ii) Compare the results of Euler's method, given in the solution to part (i) above, with 
the true solution. 


[Solution on p. 57] 


Although the results in Exercise 3 are not very accurate, the values obtained for Y, 
do follow the general pattern of the solution. It would be reassuring if this were 
always the case, but unfortunately this is not so. This can be seen in the next 
exercise, 


Exercise 4 

Apply Euler's method to find y(1) using a step length h = 0.2 for the differential equation 
dy ody 
—; + 12 + 20y =0 
de noe et Oe 


a 
with y =I and re = 1 when x = 0. Compare your results with the true solution. 
[Solution on p. 57] 


The reason for the behaviour shown by the application of Euler's method in 
Exercise 4 is discussed in the unit on numerical solutions of differential equations. 
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Summary of Section 5 
Euler's method for obtaining numerical approximations to the solution of a 
second-order differential equation is to put 
i 
dx 
and replace the second-order differential equation by a pair of first-order 


equations. Euler’s method as described in Unit 2 can then be applied to each of 
these first-order equations. 


End of section exercise 


Exercise 5 
Use Euler's method with step length h = 0.2 to find an approximation to y(1), where y 
satisfies 


@y dy 
Soe tyied 
ae ax +” 


dy 
with y = 1 and fe = 1 at x = 0. Compare your results with a graph of the analytic solution. 


[Solution on p. 58) 


6 End of unit exercises 


Exercises 1, 2 and 3 provide further practice in the main skill that this unit aims to 
teach—the solution of second-order constant-coefficient differential equations. You 
should try to make sure that you are confident on this point. If you are, you may 
like to try some of the other exercises in this section, which use the material in the 
unit, but do not concentrate solely on this central skill. 


Exercise 1 
For each of the equations below find the solution satisfying the given initial conditions. 


4» Py ody — eS 
(i) ae +9 + 5 9 with y= 2and = atx 0. 


2 
are 


ed =en 
dx? © “dx , 


i dy 
(ii) with y= 0 and = 1 at x =0. 


ay OO F dy 

(iii) ge Y= 3sin3e + 18; with y = Q, and 2 = Oate =. 
2, 

due 


[Solution on p. 58} 


di 
(iv) (u > 0); with 0= Land 2 = 2 at w= 1. 


Exercise 2 

Find the general solution of 
2 

wi) pe + 16y = 400s 4t + 8sin 4r; 


ny. doe ay eS 
(ii) ge + 4g ty = 3e 4 Bt. 
[Solution on p. 59] 

Exercise 3 


For each of the following differential equations sketch the graph of the solution satisfying 
the given initial conditions. 


i d) d 

(ye eeye0s wth lade -oersero 
xa ae dx 

. dy dy ; - dy 

(ii) wet a t Y= with y= 0 and 7 =Latx=0. 


[Solution on p. 59] 
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Exercise 4 
Suppose that x denotes the position of an object at time ¢, and that x satisfies the 
differential equation 
Px dx : 
att gy Pte te 
(i) | Sketch a graph showing the variation of position with time, a long time after the 
object is set in motion. Does the value of x ever become negative during this part of 
the motion? 
(ii) Can you make more precise the meaning of the phrase ‘a long time after the object is 
set in motion’ in part (i). 
[Solution on p. 59] 


Exercise 5 


For what values of the constant 4 can a solution other than the zero function be found for 
the boundary condition problem 


u(x) + 2ntu(x)=0; — u(0)=0, ull) =0? 
[Solution on p. 60) 
Exercise 6 
This exercise looks at an application of the theoretical ideas introduced in this unit (in 


particular, Theorem 1 of Section 3). It concerns a possible approach to the formal definition 
of the trigonometric functions. 


We define the functions S and C to be the solutions of the differential equation 
u"(x) + u(x) =0 (1) 
satisfying the conditions, 
S0)=0 S'(O)=1 
and 
c)=1 C(0)=0. 
respectively. 
(i) _ Explain why this is sufficient to define the functions S and C. 
(ii) Prove: (a) that S’ is a solution of Equation (1); (b) that S’(0) = 1 and S"(0) = 0. 


Deduce that S' = C. we ldleontinusna ded 
(iii) Use an argument similar to that in (ii) to prove that C’ = —S. aia epe na chine ee 
(iv) Find the derivative of C? + S?, and deduce that trigonometric functions in the 
C4+S=1 same sort of way, but this is 
3 sufficient to give you the 
[Solution on p. 60) flavour of the approach. 
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Appendix 


Solution to the exercise in the Introduction 


1. @ Equations (a), (c) and (d) are second-order. Equation 
(b) is first-order. 


(ii) Equations (a), (b) and (d) are linear. 
(iii) Equations (b) and (c) have constant coefficients. 
(iv) Equations (a) and (c) are homogeneous. 


Solutions to the exercises in Section 1 


1. The solutions are given in the following table: 


Auxiliary equation 


#-64+5=0 
#B+9=0 
#-4=0 
#+2M+1=0 
+4,4+8=0 
W+3=0 


=I only) (ii) 


In each case A and B are arbitrary constants. 


2. (i) The auxiliary equation is 
B+ 2wh +0? =0 
ie. +o)? =0. 


This is Case (ii), with the only root 4 = —q. The general 
solution is therefore 


y =e" "(Ax + B), 
where A and B are arbitrary constants. 
(ii) The auxiliary equation is 

5H? + 6a + Su? = 0 
thus 


—60 + /360* — 10007 
10 


=Fl-3 + 4). 


This is Case (iii) of Procedure 1.1, and the general solution is 
therefore 


3 4 a 
x= exp (— ont) (cos head + Bin; at). 
where A and B are arbitrary constants. 
3. If we can find two linearly independent solutions, then we 


can use Theorem 2 to write down the general solution since, 
with the domain (x > 0), the differential equation is of the 


y = Ae* + Be™ 

6 = Acos3t + Bsin3r 

z= Ae™ + Be-™* 

y = Ae~* + Bxe™* 
u=e~*(Acos2t + Bsin2r) 
ym A+ Bew2™2 


form specified in the Theorem. To find some solutions, we 
try y=, as suggested. fy = x* then 9 = nxt and 
2 
Sea nln Lye? Thus 
dy dy 
aby _ xo) 
cal iar 


+ y= 2x? a(n — 1)x*-? — xn 4 x 

= x"[2n(n — 1) —n +1) 

= (2n* — 3n + 1x", 
This will equal zero (for all x) if and only if 

. Int —3n+1=0. 
Teta Oa —1)(o= t= 
Hence n = $ or 1. So y = x and y = x'/? are solutions of the 
given differential equation, and these are linearly 
independent. By Theorem 2 the general solution is 


y=Ax+Bx'? — (x>0) 


4.i) The auxiliary equation is 34 — 1 — 247 = 0. 
(ii) This can be rewritten as 
27 -3+1=0 
ie. (2A-1)4-1)=0 
so A=lor}. 


(iii) By Procedure 1.1 the general solution is 
y = Ae + Bet* 
where A and B are arbitrary constants. 


5. The solutions are given in the following table: 


General solution 
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0 (only) 


In each case A and B are arbitrary constants. 


6. The general solution of the differential equation contains 
sine and cosine terms only in case (iii) of Procedure 1.1, 
when the solutions of the auxiliary equation are complex. 
The auxiliary equation is 


P+ 4h+4=0 
which has complex solutions when 


16k? — 16 <0 
is. <1 
ie, -1<k<1. 


7. No. For i 2 to apply we require that the 
coefficient ot a » never takes the value zero. But x? = 0 when 


#< t/andlia ti eetanca URW na a ede 
the domain. 


However we can say that 
y=Ax+Bx? (x >0) 

is the general solution of 
ety 


ae - 2x2 +2y=0 (x>0), 


Solutions to the exercises in Section 2 


1. (i) The associated homogeneous equation is 


p, 
Oey 0. 


Solving this gives the complementary function 

y = Acos2x + Bsin 2x, 
Gi) Trying a solution of the form y= C (C constant) in 
Sax + 4y = 8 gives 0 + 4C = 8, that is C = 2. Thus a 
particular solution is 

y=2, 


(iii) By Theorem 1 the general solution (of the 
inhomogeneous equation) is 


y = Acos 2x + Bsin2x+2. 


2. (i) For y = mx + nto be a particular solution, we require 
—2m + 2(mx + n) = 2x +3 (for all x), 


This can only be true for all values of x if m and n satisfy the 
simultaneous equations 


—2m + 2n=3 (equating constant terms) 
2m =2 (equating coefficients of x). 
Solving for m and n we obtain 
5 
m=l,n= 


z 


y=e *(Acosx + Bsinx) 
y = Ae* + Be~* 
y = Ae™ + Bxe** 


=A+Bt 


So a particular solution is 


a 
¥y z 


(ii) For y = mx + n to be a particular solution, we require 
2m + (mx + n) = 2x. 


This can only be true for all values of x if m and n satisfy the 
simultaneous equations 


2m+n=0 (equating constant terms) 
m=2 (equating coefficients of x). 


Substituting m = 2 into 2m + n =0 gives n = —4, soa 
particular solution is 


y=2e-4 


3. If we try 
y=me*, 


dy ai @y = 
thea = ~me and 7a = me '. So we need 
= 2me~* + 2me~* + me~*. 


This is satsed if 2 = Sm, that is m = = So a particular 


solution is 


y an tent 
oor al 


4. We can either try a sinusoidal solution directly or use the 
phasor method. 


Method I: We try 
y =mcos 3t + nsin 3t. 
then 


d 
& = 3ncos 3t — 3msin 31, 


dt 


— 9nsin 3t. 


So we require 
(—9m + 9n) cos 3t + (—9n — 9m)sin 3t = 


This can only be true for all values of t if m and n satisfy the 
simultaneous equations 


—9m + 9n=9 
—9m — 9n =0. 


Solving these equations gives m= —}, n = 4. So the 
particular solution is 


y = —teos3t + 4sin 31. 


9cos 3¢, 
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Method 2: We rewrite the right-hand side of the differential 
equation as Re(9e**) and try y = Re(ze*“). Using Procedure 
2.2(a) we obtain 


(-9 + 912 =9, 
i the pasticela, sclntion 
—949 -1+i 2°72 
is therefore 


y = —4cos31 + 4sin 31. 


5. Method 1: We try 
y = mcos 2t + nsin 2t. 
Then 


oY  aacon2t — 2msin 21, 
dt 


and 
d’y 
de 
‘So we need 
—8mcos 2t — 8nsin 2r = 8cos2t + 16sin 21, 


This can only be true for all values of t if m = —1 and 
n= —2. Thus a particular solution is 


= —4mcos 2t — 4nsin 2t. 


y = —cos 2 ~ 2sin 2t. 


Method 2: We write the right-hand side of the differential 
equation as Re((8 — 16i)e*") and try y = Re(ze*"), Using 
Procedure 2.2(a) we obtain 


(—4 — 4)z = (8 — 161). 
So z= —1 + 2i and the particular solution is 
y = —cos 2t — 2sin 2¢, 
6.(i) We can see that y = cos 2r is a solution of 
ey + 4y = 0, so we have an exceptional case. Two 
alternative methods of solution are given below. 
Method 1: We try 
y = t(mcos 2t + nsin2r) 
then 
& = (mcos 2t + nsin 2t) + t(2ncos 2t — 2msin2r) 
and 
dy o A 
ce (4n cos 2t — 4msin 2t) + (—4mcos 2t — 4nsin 21). 
&y. 


Substituting into ae 


+ 4y gives 4ncos 2t — 4m sin 21, so we 
require 
4ncos 2t — 4msin 2t = 2cos 2t. 


This can only be satisfied for all values of t if m = 0 and 
n= 4. So the particular solution is 


y =4tesin 20, 


Method 2: We rewrite the right-hand side of the equation as 
Re(2e?") and try 


y = Re(tze*") 
then 


¢ = Re((1 + 2it)ze%) 
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and 
= Re((—4t + 4i)ze**), 


2, 

Substituting into ay gives Re(4ize*"), so we need 
4iz = 2. 

Thus z = —4i and the particular solution is 


y = $esin 2¢, 


(ii) 2 = 1 is a solution of the auxiliary equation 

#? — 34 +2=0. So y = e* is a solution of the associated 
homogeneous equation, and we have an exceptional case. 
Following Procedure 2.3 we try 


y= mxe*. 
dy 


dy 
Then GF = ml + x)e* and 75 


= m(2 + x)e* and so 
e 
53-3 4 2y = (2-+.x)— 3 +2) + 2e)met = —me, 


Thus we have a particular solution if m = —4, The particular 
solution is therefore 


y= —4xe% 


7. In this case the general solution of the associated 
homogeneous equation is y= Ae* + Bxe*. So not only is 

y = me* going to reduce to zero if substituted, but so is 

y = mxe*. So we need to try something else to find a 
particular solution. The trick is to multiply by x again. That 
is, We try 


yamee, 

dy . @y > 
Then 7 = Gx +x yme* and a = (2x + x +2 + 2x)me*. 
So 

2, 

$Y 29 4 ym (xt 4 he +2202 + 28) + 2 me 


= 2me*. 
Thus we have a particular solution if m = }. The particular 
solution is therefore 
y=hie. 
8.3) We know from Theorem 2 that we can deal with the 
terms on the right-hand side of the differential equation 


separately, and we know from Exercise 7 that y = }x?e* is a 
particular solution for the case where the right-hand side is 
#. 


Let us now deal with the e* term. To find a particular 
solution of 


Py idy 2s 
Fe oh lar () 
we can try y = me**. This satisfies the equation if 
(4m — 4m + m)e** = e*, 
That is, if m = 1. So a particular solution of Equation (1) is 


y= 


Using Theorem 2 a particular solution of the given 
differential equation is, therefore, 


y = 4(hx7e*) — 3e** = 2x7e*— 3e3*, 
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(ii) We can find a solution by trying 
y =1+ mcos4t + nsin4t 
(or by the complex number method). 
This is a solution if (substituting in the differential equation) 


4(— 16m cos 4t — 16nsin 4t) + 9(—4msin 4t + 4ncos 4t) 
+ 100(! + mcos 4 + nsin 4t) = 9cos4r + 100. 


That is, if 
100! + (36m + 36n)cos 4t + (—36m + 36n)sin 4¢ 
= 9cos4t + 100. 
This can only be true for all values of ¢ if 


36m +36n=9 — (equating coefficients of cos) 
—36m + 36n=0 (equating coefficients of sin) 
=100 (equating constant terms). 


Solving these simultaneous equations gives | = 1, 
n=} and m=k.Soa particular solution is 


y =1+ {cos 4r + sin 41). 


9. (i) First find the complementary function, by solving 


a 
+40 = 0. 


This has auxiliary equation 
#+4=0 


which has solutions 4 = +2i. Hence the complementary 
function is 


0 = Acos2t + Bsin 2t 
where A and B are arbitrary constants, 
To find a particular solution, try 
O=mt+n. 
This satisfies the given inhomogeneous differential equation if 
A(mt +n) = 2t (for all 1); 
that is, if m=4and n=0. 
So the general solution of the given equation is 
6 = Acos2t + Bsin 2t + t. 


(ii) The associated homogeneous equation is 
ey _ dy 
dx? “dx 
Its auxiliary equation is 
347-24-1=0 


which has solutions 4 = 1 and 4 = — 
complementary function is 


Y= Ae* + Be, 
where A and B are arbitrary constants. 
To find a particular solution, try 


a -Y=0. 


4. Hence the 


Y¥=me* +n. 
This satisfies the given equation if 
3(4me**) — 2(2me**) — (me** + n) = e* + 3. 
That is, if m =4 and n= —3. 
So the general solution of the given equation is 
Y= Ae + Be* + he — 
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(iii) The auxiliary equation of the associated homogeneous 
equation is 


7 —34+2=0 
which has solutions 4 = 1 and 4 = 2. So the complementary 
function is 
y = Ae + Be™, 


where A and B are arbitrary constants. Note that in finding 
a particular solution we have an exceptional case. Both of 
the functions on the right of the given differential equation 
are solutions of the associated homogeneous equation, so for 
a particular solution we should try 


y =mxe* + nxe™, 
Then 


(1 + xymet + (1 + 2x ne 


and 


Py 


4. (2 + x)me* + (4 + 4x)ne™ 


@y 
Substituting into 34 + 2y gives —me* + ne** so we 


ae 
require 
—me* + ne™* = 2e* — Se™. 


This can only be true for all values of x if m = —2 and 
n= —5. Thus the general solution of the given differential 
equation is 


y = Ae* + Be™ — 2xe* — Sxe™. 
(iv) The auxiliary equation of the associated homogeneous 
equation is 
2 + wi + 2w* = 0, 
Thus 
A= 1-2 + \/40* — Ba) = o(-1 4). 
So the complementary function is 
y=e"(Acosat + Bsinwt). 


I will deal with the two parts of the right-hand side of the 
differential equation separately. (This is not necessary —I just 
prefer this way here.) A solution of 

= + 20 + 2u*y = w7h 
can be found by putting y = m (where m is a constant), 
which yields the particular solution 

y=th 


A solution of 


& 
£Y +208 + 20%y = Ssin 20, a) 


can be found by putting 
y =mcos wt + nsin Zot 
(or by complex numbers). Then 


oe = 2wn cos 2ot.— 2wmsin Zot 
and 
#y 
ies —4e7mcos et — 4e0?nsin oot. 
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Substituting into the left-hand side of the differential 
equation gives 

(2m + 4n)o? cos 2ot + (—4m — 2n)w? sin Zot, 
We require this to equal 5sin 2wt, which means that m and n 
must satisfy the simultaneous equations 

—2m + 4n=0 

—4m — 2n = S/a?. 
Solving these equations gives m = -4 and n= ~ sip So 
a particular solution of Equation (1) is 

1 
y= ~ Zot (2008 2oot + sin 2wt). 


Using Theorem 2, a particular solution of the given 
differential equation is 


1 C 
y= Far 2008 2wt + sin 2wt) + $h. 


Hence the general solution is 


y=e-™(Acosat + Bsinox)~ 51; Qcos2or + sin2ot) + 4h. 


10 (i) y = me**. 
(ii) y = mcos 3x + nsin3x (or y = Re(ze™*)). 


(iii) Note that the complementary function is 
y = Ae** + Be~**, so we need to try y = mxe~** for a 
particular solution. 


(iv) Note that the complementary function is 
y = Acos2x + Bsin 2x, so we must try 


y = kxsin2x + lxcos2x + mx +n 
for a particular solution (or y = Re(zxe**) + mx + n). 


11 (i) The complementary function is 
y =e *(Acosx + Bsin x). 
A particular solution is 
y=, 
So the general solution is 
y =e *(Acosx + Bsinx) +2. 
(ii) The complementary function is 
y = Ae™* + Be-™. 


For a particular solution, we try y = Re(ze*), By Procedure 
2.2(a) 
(-9-4)z= -i 


ie z= 


1 
Thus a particular solution is 
Ly 
y= 3x. 
Hence the general solution of the given differential equation 
is 
y = Ae™* + Beo** — monde 


Solutions to the exercises in Section 3 
1. () Do not be thrown by the different notation I have used 
here. The question could equally well be written; 

2 


dy . BP. dy ud 
qt 0; y=Oand 7 = 1 whent =>, 
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where y = u(t). (Letters other than y and 1 can be used for 
the variables, of course.) The auxiliary equation is 


#B+4=0. 

Therefore 4 = +2i and the general solution is 
u(t) = Acos2t + Bsin 2t. 

Now u(x/2) = 0 and so 
0=Asinx + Bcosx 


=Ax0+Bx -1. 
Therefore B = 0 and hence 
u(t) = Asin 2¢. 


Now u'(t) = 2A cos 2t and so the condition u'(n/2) = 1 gives 
1 =2Acosx = —2A. 


ie. A= —4 
‘The required particular solution is therefore 
u(t) = —4sin 20. 


(ii) The complementary function is 
y = Ae + Be* 
(see Frame 1 for details). 
A particular solution is 
y=2 
The general solution is, therefore 
y= Ae + Be* +2. 
When x = 0, y = 4 and so 4= A + B + 2. That is, 
A+B=2. (1) 


dy 20 Ca 
Now # = Act + 28e% and when x = 0, 2% = —1. Thus 


-1=A4+2B. (2) 


Solving the Simultaneous Equations (1) and (2) gives A = 5, 
B = —3. The required particular solution is therefore 


y = Se* —3e% +2. 
2. This is covered in Frame 7 of the audio-tape. 


3. Problems (j) and (iii) are initial condition problems, (ii) 
and (iv) are boundary condition problems. 


In each case, the differential equation is the same. Its general 
solution is 

u(x) = A cos 2x + Bsin 2x, 
The derivative is 

u'(x) = —2A sin 2x + 2B cos 2x. 
The solutions satisfying the given conditions are found 
below. 


(i) The condition u(0) = 1 gives A = 1. The condition 
u'(0) = 0 gives 0 = 2B. That is B = 0. The required solution 
is therefore 


u(x) = cos 2x, 

(ii) The condition u(0) = 0 gives 
0=A+0. 

The condition u(x/2) = 0 gives 
0=A-0. 


Thus B is arbitrary and A = 0. Any solution of the form 
u(x) = Bsin 2x 
satisfies the problem. 
(iii) The condition u(0) = 0 gives A = 0. The condition 
u'(0) = 0 gives B = 0. So the required solution is 
u(x) = 0. 
(iv) The condition u(—x) = 1 gives A = 1. The condition 
u(x/4) = 2 gives B = 2. So the required solution is 
u(x) = cos 2x + 2sin 2x, 


4. (i) The required solution is 

y=0 
as we have a homogeneous differential equation with initial 
conditions of the form y = 0 and 2 =0. 


(ii) The complementary function is 
y = Ae + Be-*, 

A particular solution is 
y=-8 

So the general solution is 
y = Ae + Be-* —8. 

The condition.‘y = 0 at x = 0' gives 


A+B-—8=0. (l) 
Now ose —.Be=*, and so the condition = 0 at 

dx dx 
x = 0' gives 

A-B=0. (2) 


Solving the Simultaneous Equations (1) and (2) gives 
A= B= 4, Therefore, the required particular solution is 


y=4e + 4e-* 8, 
(iii) The general solution is 

x = (At + Ble™™ 
(see Exercise 2(i) of Section 1). The condition ‘x =a at r=0° 
gives 

a=B. 


Now & = — (At + Be“ + Ae~™, and so the condition 

“dx fs 

a 7 Oat t= 0 gives 
0 = —oB + A. 


Thus A = wB = aw. The required particular solution is 
therefore 


x = alot + lew ™. 


5. (i) For first-order equations the theorem is as follows. Let 


wey + (xy = f(x) 


be any linear first-order equation for which p(x) does not 
take the value zero for any value of x. Then there is one and 
only one solution of this differential equation which also 
satisfies a condition of the form y = b when x = a. 
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We can find some solutions to the differential equation in 
(ii) by the following incautious manipulation. Separation of 
variables gives: 


fe-ft 

ee Es 

so log,y =logex+C 
ie. y= Ax. 


I say incautious because the manipulation performed is 
invalid if x or y is zero or negative. However it happens that 
the solutions found here are alright for all values of x and y, 
as I can show by checking, 


Check: If 


y=Ax 


4 
apm Ax=y 


for all real x, as required. Furthermore, if x = 0, then 
y = 0 (whatever A is). Hence all functions of the form 


y=Ax 
satisfy the given differential equation and initial condition. 


(iii) There is not a unique solution in this case. We do have a 
linear equation here. What is wrong is that the condition 
concerning p(x) is not satisfied (for p(x) = x is zero when 

x = 0). This example illustrates the fact that the condition 
p(x) # 0 in the theorem is crucial. Without it, the theorem 
does not work. 


6.i) The general solution is 
u(t) = e~*(A cost + Bsint). 

The condition u(0) = 0 gives 0 = A and hence 
u(t) = Be~* sint. 

Thus 
u'(t) = Be~™ (cost — 2sint), 


and so the condition u'(0) = 2 gives 2 = B. The required 
solution is therefore 


u(t) = 2e~* sin t. 
This is an initial condition problem. 


(ii) The general solution is 
u(t) = A cos 3t + Bsin 3t, 
The condition u(0) = 0 gives 0 = A and hence 


u(t) = Bsin 31. 
Thus 
u’(t) = 3B cos 3t 


and so the condition u’(z/3) = 1 gives 1 = 3Bcosx = —3B. 
That is B = —4, The required solution is therefore 


Oe 
u(t) = —3sin3e. 


This is a boundary condition problem. 
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(iii) We have a homogeneous equation, with initial 
conditions equal to zero. The required solution is therefore 


y=0. 
This is an initial condition problem. 


Solutions to the exercises in Section 4 


1. The phasor of cos 5t is 1 and that of sin St is —i, thus the 
required phasor is 


2=443i. 
The amplitude is 
b= /# +3 =s. 


The phase is 


Arg: = arctan (3} = 0.64. 


2. The auxiliary equation is 
B+ 2ah +07 =0 
which has roots 
A= -at/e 
The roots will be real provided «* — w? > 0. 


a. 


3. The graph is a sinusoid curve of angular frequency « (the 
amplitude and phase depend on the initial conditions). 


4. (i) We require 2* — w* to be negative and 2 to be positive 
(see table on page 36). Thus we require 


O<a<w 


(ii) (a) For this equation x = 1.5 and w = 3. Thus in this case 
0 < <«, and so the solution is a decaying oscillation. To 
ensure that the solution also fits the given initial conditions it 
must take the form shown in Figure | (starting at y = 4 with 
zero gradient). 


Initial gradient= 0 


Figure 1 


(b) For this equation x = 4 and w = 2. Thus in this case 

% > > 0 and so the solution is a combination of decreasing 
exponentials. To satisfy the given initial conditions the 
solution must take the form shown in Figure 2. 


, Initial gradient = 1 


Figure 2 


5. The solution is the sum of (/) a particular solution, which 
will be a sinusoid, and (IJ) the general solution of 


@y 4 
Gata toy=0 


Since « > 0, this general solution will be decreasing, and will 
tend to zero as x becomes large and so is a transient. Since 
« <« this transient is oscillatory. The initial conditions only 
affect the transient term in the solution. So whatever the 
initial conditions the long-term behaviour of the solution is 
the same. It is a sinusoidal oscillation, of frequency v. The 
amplitude and phase of this oscillation depend on y, a and w. 
(Example 1 on page 39 illustrates how the amplitude and 
phase can be calculated.) 
6. The general solution of the differential equation is 

y = Acos 3x + Bsin 3x. 


The angular frequency of this oscillation is 3. For the 
condition ‘y = 4 when x = 0" to be satisfied we require 


A=4, Now 4 = —3Asin3x + 3Bcos3x and so for the 


condition 2 =9 when x = 0" to be satisfied we require 3B 
= 9, that is B = 3. Thus the required particular solution is 
y = 4cos 3x + 3sin 3x. 
The phasor (see Unit 5) of this oscillation is 
z=4-3i. 
Hence the amplitude of the oscillation is 
lz] = /47 + 3? = 5, 
and the phase is 
Argz = —arctan (3) = —0.64, 
7. All the solutions of u(x) + w*u(x) = 0 are sinusoidal 
oscillations of angular frequency «. These solutions repeat 
themselves after a period = and so any solution u(x) of this 
differential equation satisfies the condition 


a 
u(x) = ulx += 
o 
for all values of x. In particular, when x = a we have 
a 
u(a) = ula + —}. 
o 


A solution can therefore never be found satisfying boundary 
conditions of the form 


u(a)=k and ua +2) 

o 

if k # . On the other hand, if k = |, then any solutions 
satisfying u(a) = k will also satisfy ula + =| =1(=k). 


Consequently the two conditions are equivalent in this case, 
and so will only place one condition on the two arbitrary 
constants in the general solution. So for k = I, the solution is 
not unique. 
8. From Example 7 of Section 1, the general solution is 

x= Ad-2+ Som 4. Bd-2- Jip, 
For this to satisfy the initial condition ‘x = a at t = 0" we 
require 

a=A+B. (1) 


Now 


= (-2 + /3pode? +9" — (2 + /3)oBel- 2 


For this to satisty the initial condition * = 0 at r= 0 we 
require 
0=(/3 —2)0A — 
Solving (1) and (2) gives 
(/3-2)4 = 2+ /3B 
=(2+4/3)(a- A). 
Thus 


1-2 fe_ (1A, 


"nee Gh 


So the required solution is 


eel ae 


x=a 


9. The solutions of the auxiliary equation 
2? + 2Paod + wo =0 
are 


j= — 2B t JAB OF = 40? 
2 
=o(-f + /p? —1). 


Hence (since f? < 1) the general solution is 
x =e" 14 cos (w,/1 — B?)t + Bsin(w./1 — B?)e 
For this to satisfy the conditions ‘x = a at t = 0° we require 
A=a 


Now 
dx pan 
5 = — Boe x 
(Acos (w./1 — B?)t + Bsin(o/1 — B?)r) 
w,/1 = pe x 
(—Asin (w,/1 — B?)t + Boos (w,/1 — B?)t). 
For this to satisty the condition S* = 0 at r= 0° we require 
«(— AB + B,/1 — f?)=0 
ie. 
AB 
Vi-F Ji-F 


So the required particular solution is 


B= 


x = ae" [cos (o,/1 — Bt + sin(o =P). 


This is a decaying oscillation, as sketched in Figure 3. 


(2+ /3)0B. Q) 
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@ Initial gradient= 0 


Figure 3 
10. To find the phasor z of the steady-state solution, we try 
y = Re(ze™), This gives 

(4(5i)? + 9(5i) + 100)z = 9 


Siz = 9. 
So 
Te 
Si 5 
Thus the amplitude of the steady-state solution is 


z= 


1 
Has 


Its angular frequency is 5. {the actual steady-state solution 


Pi 1 
is, therefore, x = feos ~ *},) 


11. If we try y = Re(ze™) we obtain the following equation 
for the phasor z of the steady-state solution: 
(=v? + 2por(iv) + w?)z = 1 
ie. 
1 
(w? — v?) + (2Bowyi° 
The amplitude of the steady-state solution is 
1 


Vio? — FF + Qfor? 


z= 


i= 


12. The nature of the solutions of the differential equation 
@ 
re — + rat +o%y=0 

depends on the nature of the roots of its auxiliary equation, 
2 + ah + w? =0. 

These are 
d= -at 2-0. 

They are real if 2? — w? is positive, complex if 2? — w? is 

negative. Let us look at the various cases in turn. 

(1) @ positive, a? — a? positive 

The roots are real, and both are negative. (The root 

—a + ,/x? — w must be negative because 2 > ,/2? — w). 
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Thus the solution of the differential equation is a linear 

combination of e** and e** with 2, and i, negative. Thus 

the typical solution in this case is a decreasing exponential. 

(ID) & negative, x? — w* positive 

Reasoning similar to that in Case (J ) shows that in this case 

the solution is a linear combination of e** and e”* with A, 

and A, positive. Thus the typical solution of the differential 

equation in this case is an increasing exponential. 

(II]) @ positive, 2? — w* negative 

In this case the roots of the auxiliary equation are complex: 
A= -ati/o?—a. 

The solutions of the differential equation are linear 

combinations of e~* cos (,/w? — a? x) and 

e~** sin (,/«? — a? x). The typical solution is decreasing and 

oscillatory. 


(IV) negative, a? — ow negative 

The formulae for the solutions are as in Case(/1J), but now 
e~™ is an increasing exponential. Thus the solutions of the 
differential equation are oscillations of increasing amplitude. 


Solutions to the exercises in Section 5 
dy 
1.@ =? 


é = —4z —Sy +cosx 
dy 
@) AF 
dz . 
ae —32sin y. 


d 
2. In each case, replace z by Ore — YeW/h, 


2 by (z,-; — 2,)/h, and x, y and z by x,, y, and z,. 
(i) For these equations we obtain 


Yes =), 


iz, — Sy, + cosx,, 
which can be rearranged to give the pair of recurrence 
relations 


Yeo =Yp + he, 
241 = —Shy, + (1 — 4h)z, + hoosx,. 


(ii) In this case we obtain 


Vers — Ve 
% %, 
and 


hai % 


= —32siny,, 


which can be rearranged to give the recurrence relations 
yp + he, 
—32hsiny,. 


Vera = 
Feet 


—02y, +z, 


The Euler approximation to y(3.2) is shown in bold type. 
2 
@) The general solution of 4% + y = 0is 


y = Acosx + Bsin x. 
The condition ‘y = 0 when x = 0’ gives A = 0. The condition 
“® 1 when x = 0" gives B = 1, So the required particular 
solution is: 

ye=sinx. 
(iii) The results are shown graphically in Figure 4. The 
general shape of the results is correct—a half cycle of an 
oscillation. However the results are not precisely accurate, 
(The values of y should not rise above 1). This inaccuracy is 
hardly surprising with so long a step-length as h = 0.2. 


Figure 4 


The Euler approximation to y(1) is shown in bold type. The 
true solution is given by 


= fille ™ — 3271), 


Figure 5 shows a comparison of this true solution with the 
Euler approximations. 


Figure 5 


The method has failed miserably to give results that are 
anywhere near reasonable, However, we need not despair 
since if we use a small enough step-length, we can get useful 
results. Clearly, we have to be careful when we use numerical 
methods. 


¥, + 0.22, 1.62, —O4y, +04 * 


12 


02 152 248 
04 2016 3.76 
2768 5610 


3.890 


The Euler approximation to (1) is shown in bold type. The 
true solution is given by 


yre*—eF +1, 


The graph indicates that the numerical results have the same 
shape as the solution curve. 


re tree eteet tee te 
+ 
: sgessetaee 
4H + ee THAR sesse 
oaspadtecattattastiogitar ctsses 
+ ; 

H 1 + 
seuss PASH ET 
sea ssaseess; : E 
reecana age aes + 
rete cat itt 


saaetesese HEE Ha 


Figure 6 
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Solutions to the exercises in Section 6 
1. @) The auxiliary equation 
#+6A4+5=0 
has roots 4 = —S, A = —1 and so the general solution is 
y = Ae~** + Be-*. 
The initial condition ‘y = 2 at x = 0’ is satisfied if 
2=A+B. (1) 
The consition = 1 at x = 0’ is satisfied if 
1=—-SA-B. (2) 
Solving Equations (1) and (2) gives A = 3 B= a The 
required particular solution is therefore 
1 
ym gllle-* — 3e-5), 
(ii) The complementary function is 
y = Ae~™* + B. 
A particular solution is y = me**, where 
4me** + 4me™* = e** 
iem= ; The general solution is therefore 


1 
at eaale Theres 


The inital conditions ‘y = 0 and “= 1 at x =0" give 


1 
O=A+B+s, 
toe ee 

~ z 


Solving these simultaneous equations we obtain A = ~% 
1 
B= ry The required particular solution is therefore 


Te see 
Tot Mane 


(iii) The complementary function is 
y = Ae™ + Be-*™, 
A particular solution of 


ee 
de 


9y = 18 


A particular solution of 
dy 7 
i 9y =3sin3¢ 

is a sinusoid with phasor z, where (putting y = Re(ze*“)) 


—18z = —3i, 
ie. Ny 
2 = si 

6 


That is 
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1 
= gains. 
The general solution is therefore 


y= Ae + BY —2— Asin 31 


4 
The conditions ‘y = 0 and 5" =0 at 1 = 0" give 
O=A+B-2 
1 
0=34-3B-5 


Solving these simultaneous equations gives A = 13, p= ii. 
‘The required particular solution is therefore 

ya hse Me~™ — 24 — 2sin 32). 
(iv) This can be solved by direct integration. We have 

dq du 1 

a-le--a*4 


We require 97 = 2 when w= 1, s0 A= 3. That is 


dé 
zo? = 
Hence 
o= {(s-2}au 
= 3u —log,u +B. 
The condition ‘@ = 1 when u = 1" gives 
1=3+B. 


‘Thus B = —2 and the required particular solution is 
0 = 3u—log,u —2. 


2. (i) The complementary function is 
y = Acos4t + Bsin4t. 
This is an exceptional case. We find a particular solution (of 


the form y = r x (a sinusoid)) by putting y = Re(zte*"). Then 
dy 
-- ze“"(1 + 4it), 

and 
> = ze*"(4i + 4i(1 + 4it)) = ze*(8i — 161). 


@ 
Set I6y = size. 


Thus the phasor z satisfies 
8iz = 4 — Bi, 


et ~ }# So a particular solution is 
a 100s 40 + 5sin 4). 
Hence the general solution is 


y = Acos4t + Bsin4t — tcos4r + jesinde 


(ii) The complementary function is 
y = Ae~™ + Bre~™, 

To find a particular solution, try 
y=lPe*+mttn 

Then 


dy 2)p-21 
G7 (Ot 2e)e-™ +m 


and 
PY 1a — 4 — e+ Ae = I= 8+ 4, 
Thus 
&y at 2 2 2) 4 422 
ettgtynk (2 — 8t + 40?) + 4(2¢ — 207) + 4¢° 
+ 4m + 4(mt +n) 
= le~™ + 4mt + 4(m +n). 

So we require 

e~* + 4mt + 4(m + n) = 3e-™ + Be. 
Therefore 2}=3, 4m=8 4(m+n)=0, That is 


jee m=) n=-2. 


? 
The general solution is therefore 
ym Aen® + Be + 3e-™ 4 24-2, 


3. (i) Since (67-4 1x1)>0 and 6>0 the graph 
has the form 


y 


Initial gradient =—2 
Figure 7 


(ii) Since (1? — 4 x 6 x 1) <0 and 1 > 0 the graph has the 
form 


Figure 8 


4.@ Since the coficient of * is positive the 
complementary function is a transient. The long-term motion 
is therefore independent of the initial conditions. The non- 
transient term will be a particular solution of the form 

x = Acos(3t+¢)+m 


where A, ¢ and m are constants. It is easiest to deal with the 
constant term m separately, For x = m to be a solution of 


dx +2 
at 
wiaaavena casi SRC 


at ene 


x=} 

z 
To find the sinusoidal solution of 
ax 
de 


we can use phasors. The phasor z of this solution is found by 
putting x = Re(ze™) (remembering that the phasor of 2sin 3r 
is —2i). So 


(i? + 31+ 4)2 = —2i 

—2i 2i(5+3i) 2 n 
ST fv gia mee 
Thus A (the amplitude) is given by 


+2 4 ae =2sin3, 


ie z= 


2 
A=(:|=—= 
is J¥ 


and ¢ (the phase) is given by 
p=Argz=n ~ arctan = 2.11 


Thus the steady-state solution is 
yee 608 (31 +211). 

2° fs 
The graph is shown in Figure 9. 


¥ 


‘The smallest value taken by x during this part of the motion 
is 


eles 
2 Ju 
This is positive, so x does not take negative values. 
(ii) ‘a long time after the object is set in motion’ means 
‘sufficiently long after the motion has started for the transient 
term to have become negligible’. We can spell this out more 


precisely if we find the transient. This is the complementary 
function, which is 


x=e [reas + asinX vi5) 


The term that determines how rapidly this transient dies 
away is clearly, e~*. Precisely when it is reasonable to 
neglect this term depends on the particular problem under 
consideration. However, it would often be eieqae to 
require that e~** < 10~*. This occurs when f > 2 log, 10* 
= 184. 
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5. The general solution of the differential equation is 
u(x) = Acos Axx + Bsin Anx. 

For u(0) = 0 we need 
O=A. 

For u(1) = 0 we then need 
0 = Bsin Ax. 


For a non-zero solution we require that B is not zero, so we 
need 


sin Ax =0. 
This occurs if A is an integer. 
6.) This is a consequence of Theorem 1 of Section 3. The 
differential equation u"(x) + u(x) = 0 is of the form specified 
in that theorem, so there is exactly one solution of this 
equation satisfying a given set of initial conditions. Hence 
there is one and only one function S satisfying the given 
definition, and so this is an adequate definition. A similar 
comment applies to the definition of C. 
(ii) (a) Since S is a solution of Equation (1): 

S'+S=0. 
Differentiating this gives 

Ss" +S'=0. 
ie, (S')"+S'=0. 
So S' is a solution of Equation (1). 


(b) Also 

s0) =1 (by the definition of S) 
and 

SOs = Aa (from Equation (1)) 


(by the definition of S). 
So the function S’ is a solution of Equation (1), that satisfies 
the initial conditions 
SO=1, (S)V)=0. 


But the function C satisfies the same initial conditions, and 
by Theorem 1 of Section 3 there is only one solution of 
Equation (1) satisfying these initial conditions. Hence 


c=S. 
(ili) C’ is a solution of Equation (1), by an argument just 
like that in part (ii), Also 

c(0) =0, (by the definition of C) 
and 

c"(0) = —C(0) 


=-l 


(from Equation (1)) 
(by the definition of C) 


Now the function —S is also a solution of Equation (1), and 
satisfies the initial conditions 


—S0)=0, -S'0)=—-1. 
Hence, by Theorem | of Section 3: 
C=-s. 


(iv) The derivative of C? + S? is 
2CC’ + 2SS' = 2C(—S) + 2SC = 0. 

Hence C? + S? must be a constant function. So 
C2 +S? = C70) + 57(0) 
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